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Introduzione 


1. Introduzione 


La configurazione del porto di Livorno e la disposizione dei varchi di accesso 
all’area portuale hanno posto alcuni ostacoli alla operatività del porto stesso. 
In particolare la mancanza di una viabilità interna al circuito doganale e, ad un 
livello superiore, allo spazio doganale, crea le condizioni per cui 1 transiti di 
merce debbano avvenire al di fuori dello spazio doganale e spesso al di fuori 
dell’ambito portuale. Tale situazione si presta ad una commistione crescente 
di correnti di traffico miste (urbano e industriale) in un zona sensibile 
dell’area portuale di Livorno, in particolare via Leonardo da Vinci. 
L'obbiettivo principale del seguente studio è quello di migliorare il livello di 
servizio degli spostamenti e dei collegamenti nella zona Nord-Ovest del porto 
di Livorno, proponendo una nuova rete viaria che permetta di suddividere il 
traffico industriale da quello urbano e valorizzando anche le zone esterne al 
porto di Livorno che a oggi non sono servite a sufficienza dalla la rete 
suburbana. 

Nel seguente lavoro, ricostruito il quadro della situazione attuale della rete 
infrastrutturale della zona Nord-Ovest del porto di Livorno, proporremo una 
nuova rete viaria di cui analizzeremo anche alcuni nodi più importanti. 

Infine verranno tratte le varie conclusioni valutando i miglioramenti ottenuti 


attraverso la realizzazione della nuova rete viaria. 
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2. Inquadramento 


2.1. Definizione della rete di interesse 


La rete attuale interessa l’area retroportuale; è delimitata a nord dal canale 
scolmatore e dalla raffineria, a ovest dal porto, a sud dal centro abitato di 
Livorno, e a est dalla periferia di Livorno. Le strade facenti parte della rete 


(evidenziata in figura 1), sono di seguito elencate: 


A - SGC FI-PI-LI 

B - Via Leonardo da Vinci 
C - Via Salvatore Orlando 
D - Via Enrico Mattei 

E - Via delle Cateratte 

F - Via Genova 

G - Via Firenze 

H - Via Federigo Enriques 
I- Via Nord 

L - Via Pian di Rota 

M - Connessione con la via Aurelia, A 12, FI-PI-LI 


I nodi di accesso/uscita della rete, facendo riferimento ai numeri riportati nella 


figura 1, sono: 


1) FI-PI-LI (Est) 

2) FI-PI-LI (Ovest) 

3) Via L. da Vinci (porto) 
5) Via Orlando 
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6) Via Mattei 

7) Via delle Cateratte 

9) Via Firenze (Sud) 

10) Via Enriques 

11) Via Nord 

12) Via Pian di Rota 

13) Connessione Aurelia 


14) Via Firenze (Nord) 


YA 


Figura 1: rete attuale 
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2.2. Analisi di traffico 
2.2.1. Rilievo dei flussi di traffico 


Per misurare i flussi di traffico è stato utilizzato in quasi tutte le postazioni un 
Radar Recorder, prodotto dall’inglese CA-Traffic Ltd. Questo strumento 
offre diversi vantaggi rispetto ai classici sistemi di misurazione, soprattutto 
per la sua facilità di utilizzo e di installazione. 

L’apparecchio è costituito da un’antenna, che emette in maniera continua 
onde elettromagnetiche ad una frequenza (24.2GHz), un mixer, e un’unità per 
la registrazione e l’elaborazione dei dati. Il segnale è inviato in maniera 
continua nella direzione in cui l’antenna è puntata. Quando viene ricevuto il 
segnale riflesso dal veicolo, questo viene comparato con un campione del 
segnale inviato. La frequenza risultante da questa operazione costituisce il 
cosiddetto spostamento (o shift) Doppler, causato dalla velocità radiale del 
veicolo verso/da il radar. Lo shift Doppler è mediato fra diversi campioni e 
processato per ottenere la velocità del veicolo. Questa viene registrata insieme 
alla lunghezza del veicolo, ottenuta misurando il tempo di occupazione del 


campo elettromagnetico (il cono rosso in figura 2). 


Figura 2: Radar Recorder 
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Il radar viene in genere messo in opera posizionandolo sui pali della 
segnaletica verticale (0 comunque su qualunque supporto disponibile) a circa 
2 metri dal suolo e con un’angolazione di 45° rispetto all’asse della strada. 

La velocità di un veicolo che si avvicina viene misurata sfruttando l’effetto 
Doppler, e poiché l’avvicinamento verso il radar non avviene secondo una 
linea radiale, il valore della velocità così calcolato viene corretto 
automaticamente con un fattore corrispondente al coseno dell’angolo di 45° 
con il quale viene posto in esercizio rispetto all’asse della strada. 

Una volta posizionato, il radar viene configurato con i dati relativi alle 
distanze della mezzeria relative alle due corsie, viene posizionata l’antenna 
con l’inclinazione indicata dal software di configurazione e viene lasciato 
registrare i dati. 

Questi dati vengono poi opportunamente tarati ed elaborati tramite il software 
Vehicle Data Analyzer 2 e Microsoft Excel. 

In alcune postazioni, a causa degli elevati volumi di traffico e la frequenza di 
code, 1 rilievi sono stati invece condotti manualmente, mediante operatore, 
solo nell’ora di punta del mattino (7:45 + 8:45). 

I veicoli rilevati sono stati suddivisi in tre classi, ovvero veicoli leggeri, 
veicoli pesanti e motocicli. Nella prima classe rientrano le autovetture e i 
piccoli furgoni dotati di asse posteriore non gemellato. Nella seconda classe 
rientrano tutti i mezzi pesanti (autotreni, autoarticolati, autobus) e i furgoni 
più pesanti (ovvero quelli dotati di asse posteriore gemellato). Nella terza 


classe rientrano tutti i veicoli a motore a due e tre ruote. 


Di seguito sono riportati l'elenco e la dislocazione delle postazioni di misura 


(figura 3): 


Rilievi mediante strumentazione RADAR 
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ID Sezione Nome Data Rilievo Intervallo orario 
01 Via Leonardo da Vinci 28/03/2017 00:00 + 00:00 
02 Via Genova 28/03/2017 00:00 + 00:00 
03 Via delle Cateratte 28/03/2017 00:00 + 00:00 
04 Viale Orlando 28/03/2017 00:00 = 00:00 
05 Connessione Aurelia 04/04/2017 00:00 + 00:00 
06 Via Firenze 04/04/2017 00:00 + 00:00 


Rilievi mediante operatore 


ID Sezione Nome Data Rilievo Intervallo orario 
07 Via Enrico Mattei 10/04/2017 07:30 + 09:00 
08 Via Firenze Nord 11/04/2017 07:30 + 09:00 
09 Via Nord 12/04/2017 07:30 + 09:00 
10 Via Enriques 12/04/2017 07:30 + 09:00 
11 Via Firenze Sud 19/04/2017 07:30 + 09:00 
12 Via Pian di Rota 20/04/2017 07:30 + 09:00 
13 a 21/04/2017 07:30 + 09:00 
14 FI-PI-LI (lato Ovest) 26/04/2017 07:30 + 09:00 
15 FI-PI-LI (lato Est) 28/04/2017 07:30 + 09:00 
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Figura 3: dislocazione delle postazioni di misura 


6 Postazione Radar 


® Postazione con Operatore 
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Di seguito sono riportati i risultati dei rilievi di traffico ottenuti mediante 


rilevazione con radar: 


Sezione 01 


Via LEONARDO DA VINCI 


Dati della Postazione 


Coordinate GPS 43°34'21.18" N 10°19'12.15" E 
Data rilievo 28-29 marzo 2017 

Periodo orario rilievo 00:00 + 00:00 

Condizioni meteo Soleggiato 


Localizzazione della Postazione 


d 


A 


Trenitalia (mi) 


©) Balzarini Moto Honda 
Costieri D Alesio S.P.A 


Immagine della Postazione 


Inquadramento 


[srecoete | SII | _0r_| 29 | 8 [srsroesi | cie | 08 | Mi | SI ] sesso 


SI'L0-00'L0 
00°L0-S#'90 
St'90-0£°90 
[c forceste | sel | de | 06 | ll |oenrsrei | ss_| 6 | & | € ]ocsosto0 
S1°90-00190 
00°90-S#"S0 
St S0-0£'S0 
0£‘SO-SI'S0 
SI'S0-00'S0 
00°S0-S#'40 
St'#0-0£ 70 
0£40-S1'40 
| OI | 5140-0070 
[_s | 9 | # [ooocsrei | 801 | e | 65 | 9 [oasi | è | è | è | 0 |oorosreo 
EG St'£0-0£'£0 
O£'£0-SI'£0 
[ci [sreroo6i | 96 | se | 0 | 8 [sroon | e | 1 | 1 | 0 [sreoooso 
[si [ooersr8i | 901 | e | #s | 6 |oonseor | _# | 1 | € | 0 |oososroo 
[uu [sesroegi | 801 | e | 9 | # |srorocor | i | 0 | | | 0 [srewworco 
[ri foesrsrgi | 001 | ee | 00 | £ |ocorsror |__| 0 | s | 0 |occosrco 
[si | «6 | ei |sreroosi | _66_| te | #0 | i [srorosor |_s_| 0 | è | 1 [sraooowo 
EC E ARI E e E EA EI 


TION HNOIZHUIC - HDUIA Gp] CIA 


Inquadramento 


<+ 
um 


UNIT 
kol Bug 
mim 


(Sl 
(| 


m 


Ci [soi | se [sreroosi | ei | e | ra | 0 [srorooor]_9_|_0 | s_|_1 |] 


TOS HNOIZAUIC - WUIA Ep "] GIA 


00°80-St'L0 


ST LO-0E LO 


ST'L0-00"L0 
00°L0-St7°90 
St#7'90-0£°90 
0£°90-S1°90 
ST'90-00°90 
00°90-St7' 50 
St 50-0E 50 
0E'S0-ST'S0 
ST'50-00°S0 
00°50-St7'70 
St'70-0E 70 
0€70-ST70 
ST'70-00°70 
00°70-S7' €0 
St E0-0E €0 
OE E0-ST'E0 
ST'E0-00°€0 
00°€0-St7'70 
St T0-0E 70 
OE TO-SI 0 
ST'T0-00°70 


10 


Inquadramento 


Sezione 02 


Via GENOVA 


Dati della Postazione 


Coordinate GPS XX N XX E 
Data rilievo 28-29 marzo 2017 

Periodo orario rilievo 00:00 + 00:00 

Condizioni meteo Soleggiato 


Localizzazione della Postazione 


(n) 


scuderia don chisciotte SON Tar Snc DI 
Zangrillo Alberto & C 


} (@ 


Ileci Di Chiama Luca ( 


Immagine della Postazione 
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Sezione 03 
Via delle CATERATTE 
Dati della Postazione 
Coordinate GPS 43°34"2.52" N 10°18'49.71" E 
Data rilievo 28-29 marzo 2017 
Periodo orario rilievo 00:00 + 00:00 
Condizioni meteo Soleggiato 


Localizzazione della Postazione 


Mariter (S.R.L.) 
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Sezione 04 
Viale ORLANDO 
Dati della Postazione 
Coordinate GPS 43°34'2.13" N 10°18'43.59" E 
Data rilievo 28-29 marzo 2017 
Periodo orario rilievo 00:00 + 00:00 
Condizioni meteo Soleggiato 


Localizzazione della Postazione 


Ekon Food 
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Sezione 06 


CONNESSIONE AURELIA 


Dati della Postazione 


Coordinate GPS 43°35'4.98 N 10°20'48.78" E 
Data rilievo 04-05 aprile 2017 

Periodo orario rilievo 00:00 + 00:00 

Condizioni meteo Soleggiato 


Localizzazione della Postazione 


Baccetti Trasporti 


©Brilli Musica 
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Sezione 06 
Via FIRENZE 
Dati della Postazione 
Coordinate GPS 43°34'47.03" N 10°20'32.81" E 
Data rilievo 04-05 aprile 2017 
Periodo orario rilievo 00:00 + 00:00 
Condizioni meteo Soleggiato 


Localizzazione della Postazione 
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00°80-St'L0 


ST LO-0E' LO 


ST'L0-00"L0 
00°L0-St7°90 
St#'90-0£°90 
0£°90-S1°90 
ST'90-00°90 
00°90-St7' 50 
St" 50-0E 50 
0£€'S0-ST'S0 
ST'50-00°S0 
00°50-St7'70 
St'70-0E 70 
0€70-ST70 
ST'70-00°70 
00°70-S7' €0 
St E0-0E €0 
OE E0-ST'E0 
ST'€0-00°€0 
00°€0-St7'70 
St T0-0E 70 
OE TO-SI 0 
ST'T0-00°70 
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[or [srecoce | ssi | ci | on | se [srsroesi | 6rI | e | o01 | _& ]srsorso 
OE LO-ST'L0 

ST'L0-00"L0 

00°L0-S#°90 

St#7'90-0£°90 

0€090-S1°90 

ST'90-00°90 

00°90-St#' 50 

St" 50-0E 50 

0 S0-ST'S0 

ST'50-00°S0 

00°50-St7'70 

St'70-0E 70 

0€70-ST70 

[srocovor |__| ci | wi | st [srerooa | e | è | di | € | srroooso 


00°70-St' 0 
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St E0-0€"€0 


OE E0-ST'€0 

[sci | se [sreroosi | 001 | oc | sor | re [sruoon] £ | © | Ss | 0 [sreoooco 

00° €0-St7'70 

St T0-0E 70 

[soi | ss [ocsrsrer] 601 | e | rai | & [ocorsror | st | I | SI | © [orcosrco 

ser [oct | rt [sreroosi | eci | 60 | sil | 0 [srorooor | st_| 1 | di | © [srevooco 


Cas HNOIZHIUIO - 2209 GIA 
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Di seguito sono riportati i risultati dei rilievi di traffico ottenuti mediante 
rilevazione con operatore manuale: 


Sezione 07 


Via ENRICO MATTEI 


Dati della Postazione 


Coordinate GPS 43°34'5.22" N 10°18'49.70" E 
Data rilievo 10 aprile 2017 

Periodo orario rilievo 07:30 + 09:00 

Condizioni meteo Soleggiato 


Localizzazione della Postazione 


Ekon Food 


Immagine della Postazione 


— 
SERE WMA (AE 20.1, 
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Via Mattei - DIREZIONE NORD 


Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 8 36 10 54 
07.45-08.00 d 47 17 61 
08.00-08.15 5 43 20 68 
08.15-08.30 di 29 16 54 
08.30-08.45 6 28 18 32 
08.45-09.00 3 39 20 62 

Via Mattei - DIREZIONE SUD 

Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 13 12 18 103 
07.45-08.00 13 99 16 128 
08.15-08.30 8 96 14 118 
08.15-08.30 8 106 20 134 
08.30-08.45 8 95 16 119 
08.45-09.00 8 103 13 124 
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Sezione 08 
Via FIRENZA NORD 
Dati della Postazione 
Coordinate GPS 43°35'16.39" N 10°20'47.71" E 
Data rilievo 11 aprile 2017 
Periodo orario rilievo 07:30 + 09:00 
Condizioni meteo Soleggiato 


Localizzazione della Postazione 


f 


v 
x 
R:2) 


© 
È 
es j 
iL 
Re, 
SS 


Immagine della Postazione 
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Via Firenze (Nord) - DIREZIONE NORD 


Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 20 297 11 328 
07.45-08.00 32 284 12 328 
08.00-08.15 20 313 14 347 
08.15-08.30 15 305 10 330 
08.30-08.45 DT DO? 24 343 
08.45-09.00 16 249 13 278 

Via Firenze (Nord) - DIREZIONE SUD 

Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 25 266 34 325 
07.45-08.00 29 257 25 311 
08.00-08.15 22, 254 21 297 
08.15-08.30 22 308 29 359 
08.30-08.45 24 284 25 935 
08.45-09.00 17 274 43 334 
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Sezione 09 
Via NORD 
Dati della Postazione 
Coordinate GPS 43°34'27.49" N 10°20'7.41" E 
Data rilievo 12 aprile 2017 
Periodo orario rilievo 07:30 = 09:00 
Condizioni meteo Soleggiato 


Localizzazione della Postazione 


lali'ASZ: 
P, 
fi O 
ale © @MBIane et Noir 
un 
da nni ARR RO 
Esso ("3 NOS Crisat Distribuzione 
<S À Elettronica SRL 
S®° 7 
D O 
e 4 
& 
<S 
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Via Nord - DIREZIONE NORD 
Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 8 46 8 62 
07.45-08.00 12 54 6 n 
08.00-08.15 i 72 6 85 
08.15-08.30 15 74 15 104 
08.30-08.45 10 79 10 93 
08.45-09.00 9 70 10 89 


Via Nord - DIREZIONE SUD 
Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 11 ST 6 74 
07.45-08.00 22 60 6 88 
08.00-08.15 10 68 8 86 
08.15-08.30 21 86 7 114 
08.30-08.45 9 80 9 98 
08.45-09.00 19 103 8 130 
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Sezione 10 
Via ENRIQUES 
Dati della Postazione 
Coordinate GPS 43°34'29.94" N 10°20'5.54" E 
Data rilievo 12 aprile 2017 
Periodo orario rilievo 07:30 = 09:00 
Condizioni meteo Soleggiato 


Localizzazione della Postazione 


@ 
< 


Crisat Distribuzione 
Elettronica SRL 


Alo'£ 
© 


[#1 EEN et Noil 
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Via Enriques - DIREZIONE SUD 


Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 1 6 6 13 
07.45-08.00 1 6 15 DI 
08.00-08.15 3 13 13 28 
08.15-08.30 3 24 18 45 
08.30-08.45 4 32 14 50 
08.45-09.00 d, 37 14 53 

Via Enriques - DIREZIONE NORD 

Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 6 13 1 20 
07.45-08.00 8 19 1 28 
08.00-08.15 6 23 8 37 
08.15-08.30 3 35 8 46 
08.30-08.45 3 29 13 41 
08.45-09.00 3 22 8 33 
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Sezione 11 
Via FIRENZE SUD 
Dati della Postazione 
Coordinate GPS 43°34'5.50" N 10°19'24.95" E 
Data rilievo 19 aprile 2017 
Periodo orario rilievo 07:30 + 09:00 
Condizioni meteo Soleggiato 


Localizzazione della Postazione 


Domenici Stefano 
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Via Firenze (Sud) - DIREZIONE NORD 


Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 S1 267 9 327 
07.45-08.00 64 290 30 384 
08.00-08.15 38 DI3 14 325 
08.15-08.30 115 390 ‘i 367 
08.30-08.45 39 266 15 320 
08.45-09.00 23 DSI 8 286 

Via Firenze (Sud) - DIREZIONE SUD 

Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 13 126 i 144 
07.45-08.00 23 172 17 212 
08.00-08.15 13 172 17 202 
08.15-08.30 16 141 14 171 
08.30-08.45 20 176 14 210 
08.45-09.00 17 148 14 179 
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Sezione 12 


Via PIAN DI ROTA 


Dati della Postazione 


Coordinate GPS 43°34'37.40" N 10°20'23.08" E 
Data rilievo 20 aprile 2017 

Periodo orario rilievo 07:30 + 09:00 

Condizioni meteo Soleggiato 


Localizzazione della Postazione 


oi 
4 
SS a È 
id 
Banca Mor 
£ Gpaschi di £ 
e O 
SO 
ZA 
£ c) 


Immagine della Postazione 
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Via Pian di Rota - DIREZIONE SUD-EST 


Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 10 96 22 128 
07.45-08.00 9 110 24 143 
08.00-08.15 4 120 24 148 
08.15-08.30 7 129 20 156 
08.30-08.45 5 109 16 130 
08.45-09.00 9 147 28 184 

Via Pian di Rota - DIREZIONE NORD-OVEST 

Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 8 86 16 110 
07.45-08.00 9 128 11 147 
08.00-08.15 6 120 19 145 
08.15-08.30 4 152 7 161 
08.30-08.45 4 125 17 146 
08.45-09.00 3 123 24 150 


37 


Inquadramento 


Sezione 13 


Via L. DA VINCI (ZONA PORTO) 


Dati della Postazione 


Coordinate GPS 43°35'9.54" N 10°19'12.16" E 
Data rilievo 21 aprile 2017 

Periodo orario rilievo 07:30 + 09:00 

Condizioni meteo Soleggiato 


Localizzazione della Postazione 


1a no @ 


Immagine della Postazione 
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Via L. da Vinci (porto) - DIREZIONE EST 
Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 1 8 13 22 
07.45-08.00 1 9 23 33 
08.00-08.15 0 3 22 a 
08.15-08.30 2 20 28 50 
08.30-08.45 0 16 33 49 
08.45-09.00 0) 28 44 74 


Via L. da Vinci (porto) - DIREZIONE OVEST 
Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 6 33 32 71 
07.45-08.00 13 54 27 94 
08.00-08.15 4 35 34 se; 
08.15-08.30 6 40 45 91 
08.30-08.45 2 29 35 62 
08.45-09.00 4 25 32 61 


39 


Inquadramento 


Sezione 14 
FI-PI-LI (LATO OVEST) 
Dati della Postazione 
Coordinate GPS 43°35'12.97" N 10°19'7.20" E 
Data rilievo 26 aprile 2017 
Periodo orario rilievo 07:30 + 09:00 
Condizioni meteo Nuvoloso 


Localizzazione della Postazione 


la no @ 
cda Vi 


Immagine della Postazione 
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FI-PI-LI (Ovest) - DIREZIONE OVEST 


Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 2 153 S1 206 
07.45-08.00 d 169 46 217 
08.00-08.15 0 171 44 215 
08.15-08.30 0 193 48 241 
08.30-08.45 1 179 47 227 
08.45-09.00 d, 165 6.1; 204 

FI-PI-LI (Ovest) - DIREZIONE EST 

Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 80 215 30 325 
07.45-08.00 46 197 44 287 
08.00-08.15 13 112 35 160 
08.15-08.30 0 102 58 160 
08.30-08.45 1 103 60 164 
08.45-09.00 5 101 53 159 
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Sezione 15 


FI-PI-LI (LATO EST) 


Dati della Postazione 
43°35'32.94" N 
28 aprile 2017 


Data rilievo 
Periodo orario rilievo 07:30 = 09:00 


Condizioni meteo 


10°20'14.21" E 


Coordinate GPS 


Nuvoloso 


Localizzazione della Postazione 


Tecnoforniture 
Industriali SRL Pr 


Floris Legnami ay [m) 
di Stenta Cindy 
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FI-PI-LI (Est) - DIREZIONE OVEST 


Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 4 12 31 107 
07.45-08.00 6 111 40 157 
08.00-08.15 3 104 35 142 
08.15-08.30 3 142 26 171 
08.30-08.45 4 128 30 162 
08.45-09.00 7 142 24 173 

FI-PI-LI (Est) - DIREZIONE EST 

Orario Due ruote Leggeri Pesanti Totali 
07.30-07.45 2 12 24 98 
07.45-08.00 2 81 33 116 
08.00-08.15 1 TO 28 104 
08.15-08.30 0 94 24 118 
08.30-08.45 0 82 39 121 
08.45-09.00 3 74 36 115 
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2.2.2. Omogeneizzazione dei dati in portate 


Inquadramento 


Non essendo disponibili per ogni punto rilevato, i dati sui veicoli a due ruote 


sono stati aggregati a quelli delle auto utilizzando il fattore di equivalenza 


appropriato (1 veicolo a due ruote = 0.5 uvp). Inoltre i dati riportati nelle 


tabelle del capitolo 2.2.1 sono stati convertiti in portate dell'ora di punta del 


mattino del giorno feriale, che nel nostro caso è 7:45-+8:45, ed è definita come 


il numero massimo di veicoli transitati in una frazione dell'ora (usualmente 4 


frazioni di 15 min, come indicato dal manuale HCM) esteso poi all'ora 


moltiplicandolo per il numero di frazioni. I risultati sono riportati nella tabella 


seguente: 
Veicoli Leggeri | Veicoli Pesanti 
(uvp/ora) (veic/ora) 
Entrante | Uscente | Entrante | Uscente 
PI FI-PI-LI (Est) 776 880 192 240 
P2 FI-PI-LI (Ovest) 396 582 160 188 
P3 Via L. da Vinci (porto) 84 242 132 180 
P4 Via Orlando 268 356 232 96 
PS Via Cateratte 292 136 12 20 
P6 Via Mattei 212 450 80 80 
P7 Via Firenze (sud) 1790 750 180 68 
P8 Via Enriques 136 156 72 SR 
P9 Via Nord 326 388 60 36 
P10 Via Pian di Rota 626 534 76 96 
P11 Connessione Aurelia 926 576 60 76 
P12 Via Firenze (nord) 1290 1316 116 96 
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2.2.3. Riequilibrio globale della rete 


Inquadramento 


Le portate ottenute devono essere ulteriormente trattate, in quanto per la 


continuità, la portata complessiva in ingresso deve essere uguale alla portata 


complessiva in uscita. Si procede quindi al calcolo delle due sommatorie: 


Di portata in ingresso = >% portata in uscita 


Portate dei veicoli leggeri entranti e Portate dei veicoli pesanti entranti 
uscenti dalla rete (uvp/ora) e uscenti dalla rete (veic/ora) 
Centroide Entranti | Uscenti Centroide Entranti | Uscenti 

FI-PI-LI FI-PI-LI 
PI (Est) 776 880 PI (Est) 192 240 
FI-PI-LI FI-PI-LI 
P2 (Ovest) 396 582 P2 (Ovest) 160 188 
Via L. da Via L. da 
P3 Vinci 84 242 P3 Vinci 132 180 
(porto) (porto) 
P4 | Via Orlando 268 356 P4 | Via Orlando 232 96 
Via Via 
BI Cateratte 200 o RO Cateratte sa 2 
P6 | Via Mattei 212 450 P6 | Via Mattei 80 80 
P7 Via Firenze 1790 750 P7 Via Firenze 180 68 
(sud) (sud) 
bal 23 136 | 156 pa 2 SH 72 52 
Enriques Enriques 
P9 | Via Nord 326 388 P9 | Via Nord 60 36 
P10 Via Pian di 626 534 P10 Via Pian di 76 96 
Rota Rota 
PII Connessione 926 576 PII Connessione 60 76 
Aurelia Aurelia 
P12 Via Firenze 1290 1316 P12 Via Firenze 116 96 
(nord) (nord) 
TOTALE 7122 6366 TOTALE 1372 1228 
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La differenza che si ha tra portate in ingresso e in uscita è pari a: 


e peri veicoli leggeri 


A= DA portata entrante — D portata uscente = 7122 — 6366 = 756 


e peri veicoli pesanti 


A= > portata entrante — ) portata uscente = 1372 — 1228 = 144 


Tale variazione, rispetto ai veicoli uscenti, assume un valore percentuale pari 


a 


e peri veicoli leggeri 


A 6 
Variazi u= ————€@< = —— =10,619 
UEAZIONE RO Y portata entrante 7122 0 
e peri veicoli pesanti 
Variazi % = a via 10,50 % 
dv Y portata entrante TAGE 


Tale valore è più che accettabile, considerando le approssimazioni che si 


hanno a livello di modellazione della rete stradale. 
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Il passo successivo è quello di riequilibrare il sistema, andando a ripartire la 


differenza di portate fra le portate uscenti e proporzionalmente ad esse. Nelle 


seguenti tabelle sono elencati 1 valori di correzione da applicare alle portate 


uscenti: 
Correzione da applicare alle Correzione da applicare alle 
portate uscenti dei veicoli leggeri portate uscenti dei veicoli pesanti 
Centroide Correzione Centroide Correzione 
PI FI-PI-LI (Est) + 105 PI FI-PI-LI (Est) +28 
P2 | FI-PI-LI (Ovest) + 69 P2 | FI-PI-LI (Ovest) +22 
Via L. da Vinci Via L. da Vinci 
Pò (borto) +29 P3 (porto) +21 
P4 Via Orlando +42 P4 Via Orlando +11 
PS Via Cateratte +16 PS Via Cateratte +2 
P6 Via Mattei +53 P6 Via Mattei +9 
P7 | ViaFirenze (sud) + 89 P7 | ViaFirenze (sud) +8 
P8 Via Enriques +19 P8 Via Enriques +6 
P9 Via Nord +46 P9 Via Nord +4 
P10| ViaPian di Rota +63 P10| Via Pian di Rota +11 
PII Connessione +68 PII Connessione do 
Aurelia Aurelia 
P12 ViaFirenze (nord) + 156 P12| Via Firenze (nord) +11 
TOTALE + 755 TOTALE + 142 
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Dopo aver applicato la correzione si ottengo i valori della tabella con le 


portate entranti e le portate uscenti equilibrate: 


Portate equilibrate veicoli leggeri Portate equilibrate veicoli pesanti 
(uvp/ora) (veic/ora) 
Centroide Entranti | Uscenti Centroide Entranti | Uscenti 
FI-PI-LI FI-PI-LI 
PI (Est) 776 985 P1 (Est) 192 268 
FI-PI-LI FI-PI-LI 
P2 (Ovest) 396 651 P2 (Ovest) 160 210 
Via L. da Via L. da 
P3 Vinci 84 271 P3 Vinci 132 201 
(porto) (porto) 
P4 |ViaOrlando| 268 398 P4 |Via Orlando 232 107 
Via Via 
PS Cateratte RE 152 PS Cateratte tz -- 
P6 Via Mattei 212 503 P6 Via Mattei 80 89 
P7 Via Firenze 1790 839 P7 Via Firenze 180 16 
(sud) (sud) 
P8 Mi 136 | 175 P8 ua 72 58 
Enriques Enriques 
P9 Via Nord 326 434 P9 Via Nord 60 40 
P10 Via Pian di 626 597 P10 Via Pian di 76 107 
Rota Rota 
PII Connessione 926 644 PII Connessione 60 85 
Aurelia Aurelia 
Pip: | Vietienze) pon. ad bip VEETIZe) ale | 40) 
(nord) (nord) 
TOTALE 7122 7121 TOTALE 1372 1370 


Le portate così equilibrate rappresentano la domanda generata ed attratta dai 


centroidi considerati nella rete. 
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2.2.4. Determinazione della matrice 0/D 


La matrice Origine/Destinazione nel nostro caso a dimensione quadrata 
12x12, pari al numero dei centroidi. Conoscendo soltanto i flussi di traffico 
entranti ed uscenti dal sistema esistono moltissime matrici O/D possibili, 
infatti manca un ulteriore vincolo che renda univoca la determinazione. 

E' solito ipotizzare una ripartizione delle portate di tipo euristico, che 
mediante una procedura iterativa ripartisce i flussi proporzionalmente e in 
modo alternato fra le portate entranti e le portate uscenti, fino a raggiungere la 
convergenza. Un'altra ipotesi è che siano nulli i flussi generati dal centroide e 
diretti verso se stesso, per cui gli elementi della diagonale della matrice 


assumeranno valori nulli in tutte le fasi del procedimento. 
Il procedimento iterativo è composta da 3 fasi fondamentali: 


1. Inizializzazione 
Avviene andando a ripartire la portata entrante dal nodo i-esimo verso 
le direzioni j, rispetto alla portata totale entrante nei nodi meno la 
portata entrante nel nodo j, proporzionalmente alla portata uscente dal 
nodo j: 
Ve,i 


A; a; =—#=<-r="moaì " 
4 i_y(Ve, i) E Ve, 


Va,j 

Con V = portata entrante V., = portata uscente 

Si ottiene una matrice di partenza, che sarà bilanciata sulle colonne cioè 
la sommatoria sulle colonne delle portate corrisponderà alla portata 


entrante nel centroide ed uscente dalla rete, mentre la sommatoria delle 
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portate sulle righe differirà dalle portate uscenti dai centroidi ed 
entranti nella rete andando a generare un differenziale A;;. 

2. Iterazione 
Nella prima iterazione il differenziale A; dovrà essere ripartito 
proporzionalmente tra gli elementi della riga generando una matrice 
bilanciata sulle righe, ma sbilanciata sulle colonne. Proseguendo le 
iterazioni continueremo a ripartire i differenziali residui, correggendo le 
righe per le iterazioni dispari e le colonne per le iterazioni pari, fino a 
raggiungere la convergenza. 

3. Test di convergenza 
AI termine di ciascuna iterazione è necessario procedere alla verifica 
del test di arresto che pone fine alle iterazioni nel caso in cui sia 
verificata la seguente disequazione: 


Max | 


Agi 
Totriga 


Aj 


, 
Totcolonna 


1< 0,01 


Nel nostro caso sono state calcolate separatamente la matrice O/D dei veicoli 


leggeri e quella dei veicoli pesanti O/D: 
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e Matrice O/D veicoli leggeri: 


Iterazione "0" 


Inquadramento 


Ae,i 
| P1 | P2 | P3 | P4 | P5 | P6 | P7 | P8 | P9 [P10|P11] P12| I 
|P1| 0 |75]30]|45|17]56 [122] 19 | 50| 71] 81]|196 -13 
|15]23|9 |[29/62|10/25|36]41]|100 16 
o0]5|2]|6]13 È 
ea [ar 26 6 o e Sw] ni 
|11|17]0 [21/46] 7 [19]27|30] 74 34 
[P6]33]21|8 |12]5]0|33]5|14|19|22] 54] 14 
| 69 [104] 40 [130] 0 | 45 |114|165|186]| 452. -34 
[P8]21]13|/5|/8|3]|10/21]0]9|12|14] 34] 15 
|P9 | 51|32|13|19] 7 |[24/51]8|0]30|34] 82 24 
61 | 24 | 36 | 14] 46|99|16|40]| 0 |65 | 158] 29 
144| 90 | 36 | 54 | 21 | 67 [146123 | 59|85| 0 | 234] 32 
[P12|200|125| 50 | 75 | 29 | 94 [203] 32 | 82 [119|134| 0 -147 
> colonna — |985 | 651 | 271 | 398 | 152 | 503 | 839 | 175 | 434 | 597 | 644 | 1472 
Aug, +|0|0{/0|0]|0/)/0|0|j0|0|0]|0|0 
NON 
Test di Arresto: 0,13 <0,01 VERIFICA 
Iterazione "1" 
z ; 
Riga AG 


PI LEZI FA PA Pe Pe Pr Pe Po PIPIL PI 


oiie-lioelielielie]ie]ie}{eEe 


(Sn) 
po 
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Au,j 2 


Iterazione "2" 


* colonna —> |985 


Au,j — | 0 


0 
0 
0 
X colonna — |988 | 656 | 273 | 401 | 153 | 506 | 839 | 176|438 | 603 | 651 | 1439 
5S|2/3]|1|3]|0|1]|4]| 6 |7 |-33 
NON 
Test di Arresto: 0,02 <0,01 VERIFICA 
DX j 
Riga AGI 
J 
1 
Ii 
(0) 
0 
0 
0 
1 
0 
0 
1 
2 
-8 
651 |271]|398|152|503|839|175]|434| 597 | 644 | 1472 
0|0|/0|0|0|0|0|0|0]|0 0 
Test di Arresto: 0,006 <0,01 VERIFICA 
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e Matrice O/D veicoli pesanti: 


Iterazione "0" 


_ 
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-22 


_ 
©) 


(0.0) 


o) 


_ 
(n 


(Sai 


Oo) 


(O) i 


_ 
N 


* colonna + 


268 


210 


201 


Nloej|pu|o|pm|wWYupr|o[|[SIN[U[ 


22, 


PAIA 


89 


N|DbUID|AI|5O|Uilr 


76 


uiwWwWlwWlolxolpailu 


58 


PB /IN[INO|N UNO | | Ul 
_ 
N 


DOD 


119 


I 
W|uwn|uuialim au 


40 | 107) 85 | 107 


Au, > 


Iterazione "1" 


Test di Arresto: 


0,13 <0,01 
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P10]1 17/13/12/7)1|5/5]|3j/2]|0%)|5]| 6 76 0 
PIL 13° | 10 [190 | 590 | | del 35] 2. 96 #06] 25 60 0 
P12f 26) 2018/11] 2|]8|7]|5]|3|9)|7]|0 116 0 
X colonna + |264|210|201|108| 22 | 90 | 77 | 59 | 40 [108] 86 (108 
Au, -—>-|4|0|0|1|0]|1 1 1|O]|1 1 1 
NON 
Test di Arresto: 0,015 <0,01 VERIFICA 
Iterazione "2" 
x ; 
Riga Aci 
I i 
3 7 191 -1 
3 5) 160 0 
2 4 132 0 
4 7 232 0 
0 1|0|0 12 0 
I > 80 0 
3 0|8 | 5 180 0 
1 4.2 72 0 
Ì Sil 20 60 0 
1 5 ed 02 76 0 
il 4|3)|2 60 0 
2 CA 116 0 
» colonna — |268|210|201|107| 22 | 89 | 76 | 58 | 40 [107] 85 | 107 
Au, —-|0|0|0|0|0]|0]|0|/0|0|0}|;0]|0 
Test di Arresto: 0,005 <0,01 VERIFICA 
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2.3. Modello della rete attuale 
2.3.1 Il software AimsunIM 


Il programma si può considerare una piattaforma integrata di modelli atti a 
replicare 1 fenomeni connessi alla circolazione dei veicoli, 1 suoi algoritmi 
permettono di modellare il comportamento su scala microscopica e offrire 
così una rappresentazione verosimile del reale comportamento degli utenti in 
strada. L’analisi può andare dalla simulazione della singola intersezione fino 
alla simulazione di intere reti urbane, permettendo l’analisi microscopica, 
mesoscopica e macroscopica. Questo programma si interfaccia con GIS, 


CAD, altri simulatori e software di pianificazione. 


2.3.1.1. La modellazione dello scenario 


Lo scenario di simulazione viene realizzato sagomando la rete o l’intersezione 
su un disegno di fondo. Lo scenario è composto da elementi a senso unico 
chiamati section, che costituiscono gli archi viari, i quali si collegano tra loro 
in incroci che prendono il nome di node. Una strada a due corsie per senso di 
marcia verrà quindi riprodotta usando due section, una orientata in senso 
opposto all’altra. L'immagine di fondo su cui verrà modellata la rete o 
l’intersezione può essere in formato vettoriale DWG o DXF, ma anche in 
formato JPEG. Strumenti appositi permettono di assegnare le opportune unità 
di misura e di scalare l’immagine o il disegno. Una volta modellate, le section 
vengono unite in un node che può essere un incrocio non semaforizzato, 
semaforizzato o a rotatoria. La disciplina dell’incrocio (regole di precedenza, 


ciclo semaforico, manovre permesse) vengono stabilite dall’utente. Quando il 
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modello è realizzato si impostano le condizioni del traffico, le quali possono 


essere inserite nel programma in due modi: 


e mediante la determinazione di uno stato di traffico (Traffic State) 


e mediante l’assegnazione di una matrice Origine/Destinazione (Traffic 


Demand) 


Con il primo metodo si conosce il flusso del traffico in ingresso nel sistema e 
ogni percentuale di svolta: è noto quindi ad ogni incrocio quanti veicoli in 
percentuale proseguono a dritto, quanti svoltano a destra e quanti svoltano a 
sinistra. Questo metodo si ottiene come risultato di un’assegnazione 
precedente, o da una solida e capillare campagna di rilevamenti, o infine come 
ipotesi sperimentale imposta dall’utente. I veicoli vengono generati nei punti 
di ingresso della rete in maniera casuale e altrettanto casualmente si 
distribuiscono per la rete, rispettando però l’assegnazione dei valori di svolta 
stabilita dall’utente. I veicoli non conosco l’intero percorso ma solo il 
prossimo movimento di svolta. 

Il secondo metodo prevede l’assegnazione di una matrice O/D. L’utente deve 
creare i centroidi, vale a dire 1 punti da cui si originano i veicoli e verso cui i 
veicoli sono diretti. Quindi viene impostata la matrice decidendo per un 
intervallo di tempo desiderato quanti veicoli generati da un determinato 
centroide si dirigono verso uno specifico centroide di destinazione. Questo 
metodo è stato impiegato nel corso della ricerca oggetto del presente 


elaborato. 
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2.3.1.2. La modellazione delle manovre dei veicoli 


I veicoli nella rete seguono i modelli comportamentali del Car Following, del 
Lane Changing, e del Gap Acceptance. Le vetture subiscono condizionamenti 
nel loro comportamento di guida in base all’ambiente circostante, alla 
segnaletica, alla geometria della strada e alla presenza e al comportamento dei 
veicoli vicini. 
Il programma suddivide la singola simulazione in intervalli di tempo, definiti 
cicli di calcolo o Simulation Step, la cui ampiezza è compresa tra 0,1 s e 1,5 s. 
AI termine di ogni ciclo di calcolo il programma calcola per ogni veicolo 
presente nella rete la posizione e la velocità e le aggiorna in base al seguente 
algoritmo: 

If (it is necessary to change lanes) then 

Apply Lane-Changing Model 

Endif 

If (the vehicle has not changed lanes) then 

Apply Car-Following Model 

Endif 


Al termine dell’aggiornamento vengono introdotti nella rete i veicoli previsti 
per il successivo ciclo di simulazione e si calcolano i tempi di arrivo per i 


veicoli del turno successivo. 
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Vengono brevemente illustrati i modelli impiegati: 


Car Following: questo modello stabilisce la velocità longitudinale del 
veicolo in base alla posizione e alla velocità degli altri veicoli, alla 
velocità desiderata e al tipo di guidatore. Questo modello possiede due 
componenti: accelerazione e decelerazione. La prima rappresenta 
l’intenzione di un veicolo di raggiungere una certa velocità desiderata, 
mentre la seconda riproduce le limitazioni imposte dal veicolo che lo 


precede quando il primo cerca di guidare alla velocità desiderata. 


Lane Changing: questo modello stabilisce se per un veicolo è possibile 
e conveniente effettuare un cambio di corsia. La decisione di cambiare 
corsia dipende dalla decelerazione che la manovra induce nei veicoli a 
monte, dalla distanza dal termine della section, dal confronto tra il 


traffico nella corsia corrente e quello nella corsia desiderata. 


Gap Acceptance: il modello valuta la possibilità per un veicolo di 
effettuare una manovra in un contesto in cui deve dare precedenza ad 
un’altra corrente veicolare. Si determina se un veicolo in approccio ad 
una intersezione da una strada secondaria può o non può attraversarla in 
funzione dei parametri dei veicoli sulla strada principale (posizione e 
velocità). Questo modello prende in considerazione la distanza del 
veicolo dall'ipotetico punto di collisione, la velocità e la sua 
accelerazione. In seguito determina il tempo necessario al veicolo per 
prendere una decisione ed attraversare l'intersezione che è anche 
funzione del livello di rischio caratteristico di ciascun guidatore. I 


parametri che hanno una maggiore influenza in questo modello sono 
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l’accelerazione del veicolo (a cui è legato il tempo di manovra), il 
massimo tempo di attesa che il conducente è disposto ad aspettare 
prima di tentare un’immissione azzardata (give way time), e il tempo di 
reazione (reaction time), che influenza il margine di sicurezza, cioè 
quel lasso di tempo aggiuntivo che il guidatore somma al tempo stimato 
per la manovra in modo da essere certo di riuscire ad oltrepassare in 


sicurezza il punto di conflitto. 


2.3.1.3. Le classi veicolari 


Aimsun®M dispone di alcune classi veicolari predefinite (auto, mezzi pesanti, 
autobus, pedoni) ed è possibile crearne di nuovi con le caratteristiche 
desiderate. Tutte le tipologie veicolari si possono definire in termini di: 

e caratteristiche geometriche 

e visualizzazione 2D e 3D per la restituzione grafica 

e parametri cinematici e dinamici (accelerazione, decelerazione, 

frenata di emergenza e tempo di reazione) 

In un approccio di simulazione dinamica, le manovre dei veicoli sono 
modellate, come visto precedentemente, utilizzando il modello di ‘“car- 
following”, ‘“lane-changing” e ‘gap acceptance”’. Questi modelli di 
comportamento del veicolo sono in funzione dei vari parametri che 
consentono la modellazione dei vari tipi di veicoli. L’utente può impostare 
questi parametri a seconda delle caratteristiche del traffico da riprodurre, del 


paese e dell'ambiente che può essere rurale o urbano. 
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2.3.1.4. L'elaborazione delle simulazioni 


Modellata la rete o l’intersezione e impostati i parametri dei veicoli, viene 
creato ciò che il programma definisce Scenario, vale a dire un insieme di più 
simulazioni della stessa situazione. In questa fase il programma impone di 
scegliere: la matrice O/D (o le matrici O/D se si utilizzano matrici per veicoli 
diversi o per fasce orarie diverse), un valore per il simulation step, il tempo di 
reazione davanti ad uno stop, il tempo di reazione davanti ad un semaforo e 
un tempo di warm up (che è il lasso di tempo dall’inizio della simulazione di 
cui il programma non tiene conto per la valutazione dei tempi di attesa, della 
lunghezza delle code, dei tempi di viaggio e di ogni altro dato di uscita) per 
dare modo alla rete di caricarsi adeguatamente. 

In questa parte viene anche imposto un valore di look ahead, vale a dire il 
numero di svolte successive che un veicolo deve fare. Normalmente ogni 
veicolo conosce solamente la manovra da compiere al termine del ciclo di 
calcolo, ma imponendo look ahead più alti (da 2 in su) si permette al veicolo 
di conoscere anche una o più manovre successive alla prossima che deve 
effettuare. 

Quando tutti questi parametri sono impostati viene creato un esperimento 
(Experiment) composto da simulazioni (Run) il cui numero viene deciso 
dall’utente, ciascuna delle quali carica la rete in maniera casuale e diversa 
dalle altre run. La disciplina degli arrivi può essere scelta tra le varie opzioni 
disponibili: esponenziale, normale, uniforme, A.S.A.P. Solitamente si opta 
per la distribuzione degli arrivi esponenziale. 

Mediante la funzione Average (media) i valori dei singoli dati in uscita 
ricavati dalle varie run vengono espressi in un unico valore globale, che 


costituisce il dato finale ricercato tramite la simulazione. 
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2.3.1.5. La restituzione dei risultati 


Ogni simulazione può essere eseguita in tempo reale, mostrando la rete con le 
vetture che vengono generate e si muovono attraverso di essa. La restituzione 
grafica può essere bidimensionale o tridimensionale. I risultati restituiti da 
Aimsun®M sono molteplici tra cui: il tempo di viaggio, il tempo di attesa, la 
lunghezza massima e media delle code, l'emissioni di anidride carbonica, 
l'entità del flusso, l'indice di saturazione, la densità veicolare. Questi risultati 
sono riportati sia per l’intera rete, sia per la singola section, per ogni 
simulazione o globalmente per l’intero esperimento. Possono essere esposti 
graficamente sul modello stesso mediante la sovrimpressione di testo o la 
colorazione. Alla rappresentazione grafica si affianca la raccolta dei dati in 
grafici o tabelle. Il programma permette di ricavare indicatori statistici dai dati 
ottenuti, come media e deviazione standard, o di confrontare i dati provenienti 
da più simulazioni rappresentandole sullo stesso grafico o appaiando 1 risultati 


nelle tabelle. 
2.3.1.6. Calibrazione 


Il programma Aimsun!M permette di creare modelli generali che devono 
essere adattati alla situazione esaminata mediante un’opportuna calibrazione. 
Il processo di calibrazione avviene seguendo un preciso schema, (riportato in 
figura 4): si sceglie un parametro legato alla situazione reale o di progetto (ad 
esempio il tempo di viaggio da un capo della rete ad un altro, o come nel caso 
del presente elaborato le portate di alcuni punti interni alla rete), e lo si 
confronta con lo stesso parametro dato come risultato da AimsunTM. Se la 
differenza tra i parametri reali e quelli ottenuti attraverso il programma 


soddisfa la verifica degli indici di discrepanza il modello si ritiene calibrato. 
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Dati iniziali 


(volumi di traffico, geometria) 


Modello di 
microsimulazione 
Valori in uscita 
simulati 


Modello di calcolo 


analitico o empirico 


Valori in uscita 
calcolati 


Verifica degli 
indici di discrepanza 


Modifica dei soddisfatta? 


parametri 


Calibrazione accettata 


Figura 4: schema di procedura della calibrazione 


Esistono vari indici di discrepanza per valutare la bontà della calibrazione: 


1. Il più noto e largamente impiegato indice di discrepanza tra modello e 
situazione reale è la Radice dell’Errore Quadratico Medio (Root Mean 
Squared Error- RMSE), intesa come radice quadrata dell’Errore 
Quadratico Medio (Mean Squared Error- MSE). Definito x; il valore del 


parametro stimato di un elemento 1 e y; il valore omologo osservato dal 
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modello per lo stesso parametro riferito allo stesso elemento, il valore 


di RMSE è dato dalla seguente espressione: 


Con n pari al numero degli elementi. 

Si passa a valori assoluti operando una normalizzazione dell’argomento 
del radicando, vale a dire dividendo la differenza tra i parametri 
confrontati per la media del parametro osservato Y. Si ottiene in questo 
modo la Radice dell’Errore Quadratico Medio Percentuale (Root Mean 


Squared Percent Error- RMSPE), la cui espressione è la seguente: 


PSE 


Di 


RMSPE = 


Il risultato dell’RMSPE è un numero puro ed esprime una stima 


dell’errore totale percentuale. 


. Altro indice diffuso è il coefficiente di correlazione lineare di Pearson r 
che indica quanto la relazione tra valori osservati e omologhi valori 


stimati si avvicina alla linearità. 


1 x (x: - DO: -7 


TYT= —_— rr Ye] 
1 
° A [n.@-2? L.01-° 


In cui ai simboli già noti si aggiungono X e y, intesi come valori medi 


dei parametri stimato x e di quelli osservato y. 
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3. Theil ha sviluppato una ulteriore misura della bontà di adattamento 
chiamata "Coefficiente di disuguaglianza di Theil", che risulta più 
sensibile ed accurato di quelli elencati precedentemente. 

Il coefficiente di discrepanza di Theil (U) ha la seguente espressione: 


1 2 
n ji - x) 


[1 g si È 


U = 


Perché 11 modello possa ritenersi calibrato viene proposto per ciascun indice 
un valore di soglia che segna lo spartiacque convenzionale tra un buon 


adattamento e un cattivo adattamento dei dati. 
RMSP < 0.20 (o meglio 0.15) 
r > 0.8 (correlazione lineare forte) 


U < 0.3 (discrepanza trascurabile) 


Nel caso in cui al termine della prima simulazione, ottenuti quindi i dati del 
programma e messili a confronto con i dati equivalenti osservati o calcolati, le 


condizioni di calibrazione non fossero rispettate, è necessario modificare i 


parametri della simulazione, elencati di seguito dal più al meno influente: 


1. Simulation Step: Indubbiamente il parametro dalla maggiore influenza 
e quindi il primo da modificare per cercare una soluzione di 
calibrazione accettabile. Più basso è il valore imposto, più gli utenti 


prenderanno decisioni rapidamente e avranno maggiore prontezza nelle 
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manovre, per cui sapranno sfruttare molti più varchi nelle correnti di 
disturbo e questo comporta tempi di attesa minori. Di contro un elevato 
simulation step, coincidente con un elevato tempo di reazione, 
porterebbe a tempi di attesa più ampi, o addirittura alla congestione di 
un ramo per incapacità degli utenti di immettersi. 

. Reaction Time at Stop: Il tempo di reazione allo stop è il secondo 
parametro per importanza. Esso è il tempo necessario perché un veicolo 
fermo davanti alla linea di stop (o nel caso delle rotatorie alla linea del 
dare precedenza) riconosca un varco nella corrente di disturbo e lo 
sfrutti per immettersi. Come è intuibile, tempi di reazione allo stop 
bassi indicano utenti molto reattivi e questo porta a ridurre i tempi di 
attesa, mentre tempi di reazione più alti indicano utenti scarsamente 
reattivi con conseguente aumento dei tempi di attesa e lunghezza delle 
code. 

. Look ahead: Questo parametro indica il numero di manovre 
consecutive che l’utente sa di dover fare. Di base ogni veicolo conosce 
solo la decisione che prenderà al termine del ciclo di simulazione, 
quindi sa solamente cosa farà alla svolta successiva. Aumentando il 
valore del look ahead il veicolo conosce le svolte e le manovre che 
dovrà compiere non solo al termine dell’attuale ciclo di calcolo, ma 
anche del successivo o dei successivi. Questo gli permette di prepararsi 
alle manovre con anticipo riducendo i tempi di attesa. Di contro 
diminuire il valore del look ahead o disattivarlo porta a tempi di attesa 
più lunghi. 

. Parametri dei veicoli: Il quarto parametro da cambiare nel caso i 
risultati fossero ancora lontani dall’offrire un modello calibrato sono i 


parametri dei veicoli, in particolar modo la velocità desiderata, la 
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massima accelerazione, la decelerazione ordinaria e la massima 
decelerazione a influire sui tempi di attesa. 

5. Queuing up speed e queuing leaving speed: Gli ultimi due parametri, 
spesso accoppiati e considerati come un solo parametro, riguardano 
prettamente le code ed esprimono rispettivamente la velocità sotto la 
quale un veicolo è considerato in coda (Queuing up speed) e la velocità 
oltre la quale si ritiene che il veicolo abbia lasciato la coda. Modificare 
le velocità che definiscono la coda ha ripercussioni di lieve entità sui 


tempi di attesa. 


2.3.2 Simulazione con Aimsun®M 


In questa fase dello studio si è creato il modello relativo allo stato attuale per 
poter effettuare la simulazione. 

La rete con le sue caratteristiche geometriche viene modellata su Aimsun!M 
basandosi sulle ortofoto e sulla Carta Tecnica Regionale messe a disposizione 
dalla Cartoteca della regione Toscana. Si è disegnato 1 nodi e i rispettivi archi 
della rete assegnando ad ognuno di essi le caratteristiche geometriche e 
tipologiche alle quali è associata una determinata velocità di percorrenza degli 
archi e la capacità. In seguito si creano le classi veicolari e si inseriscono i 


centroidi. Un'immagine del modello è visibile nella figura 5: 
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QUADRO 1 


QUADRO 4 


QUADRO 3 


QUADRO 2 


Figura 5: modello della rete e sua suddivisione 


A questo punto le matrici trovate nel paragrafo vengono inserite in AimsunTM 
che, utilizzando l’algoritmo di Frank e Wolfe, ha determinato la distribuzione 
dei flussi negli archi della rete in condizioni di equilibrio deterministico, sia 
complessivamente (in termini di unità veicoli passeggeri - uvp) sia per i soli 
veicoli leggeri e i soli veicoli pesanti. Inoltre ha aggiornato le matrici O/D in 


modo che lo scarto tra 1 flussi stimati dal modello e quelli misurati fosse 
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minimo. Di seguito si riportano le matrici O/D aggiornate dei veicoli leggeri e 


pesanti e la matrice O/D dei veicoli complessivi espressi in uvp (1 veicolo 


leggero = 1 uvp; per uniformare i veicoli pesanti, si è utilizzato 1 veicolo 


pesante = 2 uvp). 


Matrice veicoli leggeri 


Centroide | P1 | P2 | P3 | P4 |P5| P6 | P7 | P8 | P9 |P10|P11| P12 
PI O |312| 88 |217|21| 85 | 52| 8 |16|22| 8 | 33 | 860 
PZ 288| 0 |31|78|7|30|17|]3|5]|7]|2]| 1 | 480 
P3 88|34|0|22/3/9|5/1/1)2|]|1|3 [MS 
P4 81|30|8 | 0 |16|65|41| 6 |12 17] 6 | 25 | 307 
pa 92 |34| 9 |13|]0|74|46| 7 |14|19| 7 | 29 | 344 
P6 84 |32|9 |10|17| 0 |42| 6 |12|18| 6 | 27 | 261 
P7 101| 36 | 9 | 94 [21] 90 | 0 |126|251|355|143| 564 | 1788 
P8 2|0]0|6|0|1]|14|0]|4]|6 | 2 | 10 [246 
P9 9|3]|1]|16|2|8|84|14| 0 |37|14]| 58 | 245 
P10 23|8 | 2|31|5]|20]|201|31|63| 0 |35 | 14 | 560 
P11 2|0|46|1|6]|72|11|22|31]| 0 | 597 | 793 
DE 23 | 8 | 2|75|5]|24]|272|42|84|118|532| 0 |1186 

796 | 498 | 158 | 608 |97 | 411 | 847 | 253 | 485 | 631] 755 | 1497 | 7038 

Matrice veicoli pesanti 
Centroide [MPI] P2 P4 P6|P7|P8|P9| P10 | P11]|P12 

0 |78]70/37|/1)7/2|/1j12]|1 1/55 

P2 8]0|502z/1|5/2|/1]1]1|]0)1 MO 
65 |4| 0) 22|/1]41]1]1]11|]0 1] 

P4 48 | 33 | 29) 0 |9]|48/17|10]7]| 16 | 4 | 9 | 230 
1|1|1])1]0|2|1]0|0o]1|0])0jM 

P6 20 |14|/12) 6 |4|0|6|/4|3]7]2]|4 MS 

e ia {3 [4 16/0 [29/22] 46|/13)25 MU 

P8 1|1|1]|6|0|1|13|0|6|12j4]|7 |M 
1|1|1])5]/0|1|]95|0]9|2]|5/5553 1 

P10 2|1|1]|6]|0|2|15|/9|7]0|4]|° MS 
1/1|/0])5/0|1|714|/4)7]/0]/31 MU 

P12 2|1]1]9]|1/3]|19[11|9]/19]|4]0 MW 
177 139 80 |92 | 76 120 MM 01] 1327 
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Matrice veicoli complessivi 


Flussi di veicoli leggeri [veic/h] 
Quadro 1 


Via Leonardo da Vinci 
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Flussi di veicoli pesanti [veic/h] 
Quadro 1 


Via Leonardo da Vinci 


Via Leonardo da Vinci 


Via Genova 


46 


Via Orlando 
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Flussi di veicoli pesanti [veic/h] 
Quadro 2 


Via Leonardo da Vinci 


Via Genova 


45 


Via Orlando 


Via Enriques 


Via Nord 


Via Firenze 
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Flussi di veicoli pesanti [veic/h] 
Quadro 3 


Via Pian di Rota 


n — 


Flussi di veicoli leggeri [veic/h] 
Quadro 4 


Connessione con Aurelia, 
A12, FI-PI-LI 


Via Firenze 0 
0 
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Flussi di veicoli pesanti [veic/h] 


Quadro 4 


Men 


Via Firenze 


"= 
(ea) 


Connessione con Aurelia, 
A12, FI-PI-LI 


Via Firenze 


2.3.3. Calibrazione e validazione del modello 


La calibrazione del modello è stata ottenuta confrontando le misure di 


traffico, rilevate e stimate, di alcuni archi interni alla rete: 


e ViaL. da Vinci 
e Via Genova 
e ViaFirenze 


Con lo stesso metodo visto nel paragrafo 2.2.2. sono state calcolate le portate 


di questi punti che si riportano nella tabella seguente: 
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Veicoli Leggeri Veicoli Pesanti 
Sud Nord Sud Nord 
Via L. da Vinci 712 566 128 172 
Est Ovest Est Ovest 
Via Genova 672 280 132 20 
Sud Nord Sud Nord 
Via Firenze 648 1204 108 60 


Sono stati necessari numerosi tentativi prima di riuscire a trovare i valori dei 


parametri del programma che permettessero una calibrazione soddisfacente. 


Nella figura 6 sono riportati i soli parametri modificati rispetto allo stato di 


default: 


Main | Behaviour ReactionTime | Arrivals | RouteChoice | Veriables | Polides | Attributes | 


Simulation Step (in seconds) 


Esse d _ ___—_—_—__________—_—_—_—_—_—_—_— 


‘Reaction Time (Related to Simulation Step) 


(Fixed (Same as Simulation Step) 


Reaction Time at Stop (in seconds) — — 


C Variable (a Multiple of Simulation Step) 


Reaction Time at Stop: |0.5 


Reaction Time at Traffic Light: [0.5 


Figura 6: parametri dei tempi di reazione 
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Inquadramento 


Si esegue quindi l’esperimento composto da quattordici simulazioni, e quando 
ogni simulazione è stata eseguita si richiede al programma di calcolare 
l’average, vale a dire di ricavare i valori globali dell’esperimento, che sono 


riportati nella tabella seguente: 


Valori Stimati 


Veicoli Leggeri Veicoli Pesanti 
Sud Nord Sud Nord 

Via L. da Vinci 655 613 121 183 
Est Ovest Est Ovest 

Via Genova 550 440 113 43 
Sud Nord Sud Nord 

Via Firenze 734 1087 92 88 


I risultati ricavati si pongono a confronto con i risultati rilevati per 
determinare gli indici di discrepanza e validare il modello. I risultati del 


confronto sono riportati nella tabella seguente: 
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Veicoli leggeri Veicoli pesanti 
Y X (da Y X (da 
(osservato) | Aimsun®M) | (osservato) | AimsunTM) 
3 Su Sud 712 655 128 121 
a E 566 613 172 183 
Via Genova Est 672 550 132 113 
Ovest 280 440 20 43 
ui Sud 648 734 108 92 
Via Firenze 
Nord 1204 1087 60 88 
RMSPE 
r 
U 


Si può vedere che 1 valori limite sono tutti rispettati, in misura più o meno 
segnata. Pertanto si può ritenere che il modello di Aimsun®M sia calibrato e 
validato. 
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3. Progetto di una nuova rete viaria a Nord-Ovest del porto di 


Livorno 


3.1 La nuova rete 


La viabilità del nuovo sistema è stata realizzata con lo scopo di : 


e Creare due direttrici di traffico, una per il traffico industriale diretto 
al porto e una per il traffico urbano diretto a Tirrenia. 

e Valorizzare le attività che hanno la loro sede nella zona Nord-Ovest 
del porto di Livorno migliorando i collegamenti. 

e Spostare il limite portuale, che ad oggi coincide con via L. da Vinci, 


ampliando la zona retroportuale. 


Nella figura seguente è rappresentato il modello della rete di progetto del 


nuovo sistema: 
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-————— Progetto "Raccordo" 


| -—— Strade attuali di connessione 


sez Viabilità di Progetto nuovo sistema 
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Per poter raggiungere gli scopi sopra elencati sono state ipotizzate alcune 


modifiche alla rete attuale: 


e Nel nodo G è stata inserita una rotatoria che divide la viabilità ora 
esistente destinando le strade di via Orlando e di via L. da Vinci al 
traffico industriale e le strade di via Mattei e di via Genova al traffico 
urbano. 

e Destinando il ponte di via Genova al solo traffico urbano è stato 
necessario ipotizzare, per il traffico industriale, una nuova direttrice che 
collegasse il nodo G a via Firenze. Questo è stato ottenuto con 
l'inserimento di un arco viario dal nodo G al nodo H, che permette di 
collegarsi a via Enriques attraverso un nuovo ponte parallelo a quello 
già esistente in via Genova. 

e Inoltre con l'inserimento dell'arco G-H si è potuto realizzare una nuova 
rete suburbana che valorizza le zone prima poco servite. 

e Queste modifiche apportate alla rete fanno sorgere il problema di 
connettere sia il traffico urbano che dalla città di Livorno si dirige verso 
Tirrenia, sia il traffico industriale proveniente dalla variante Aurelia e 
dalla raffineria diretto verso il porto. Il problema viene risolto con il 
prolungamento di via Enriques che dal nodo H, dopo aver scavalcato la 
ferrovia con un ponte, da una parte colleghi con un nuovo arco viario il 
traffico urbano a Tirrenia e dall'altra colleghi attraverso una rotatoria 
aerea nel nodo I il traffico industriale alla FI-PI-LI e alla zona portuale. 

e Peralleggerire una buona parte di via Firenze dal traffico industriale e 
valorizzare le zone adiacenti alla raffineria adesso abbandonate, è 


prevista la realizzazione di un nuovo arco viario che colleghi la 
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rotatoria già esistente nel nodo A al nodo H mentre nel nodo B, tra via 
Firenze e via Enriques, è previsto l'inserimento di una rotatoria al posto 


dell'intersezione semaforica. 
3.1.1. Modello del nuovo scenario 


Partendo dal modello dello stato attuale, in base alle modifiche espresse nel 


paragrafo precedente, è stata creata su Aimsun®M la nuova rete (figura 7): 


Ù 50) 
. PI 
è ° 
o A 
Pg ® 5 a Pi? 
PO (1) Ln ° PI 
°° 
P2 è Pes [ 
ed 
P3 . 


® 
P10 


% 


® ° 
p4 0 o) 
P6 ps P7 


Figura 7: modello della rete di progetto 
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3.1.2. Ridistribuzione flussi di traffico 


Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


Per effettuare la ridistribuzione dei flussi veicolari nella nuova rete è stato 


necessario partire dalle matrici O/D utilizzate nella simulazione su AimsunTM 


della rete attuale presenti nel paragrafo 2.3.2.. 


Da queste matrici si ottengono 1 flussi di traffico entranti e uscenti dello stato 


attuale, riportati di seguito: 


Flussi di traffico entranti e uscenti Flussi di traffico entranti e uscenti 
della rete stato attuale dei veicoli della rete stato attuale dei veicoli 
leggeri (uvp/ora) pesanti (veic/ora) 
Centroide Entranti | Uscenti Centroide Entranti | Uscenti 
P1 | FI-PI-LI (Est) 862 796 P1 | FI-PI-LI (Est) 201 228 
FI-PI-LI FI-PI-LI 
P2 (Ovest) 479 499 P2 (Ovesi) 171 179 
P3 Via L. da Vinci 169 159 P3 Via L. da Vinci 142 169 
(porto) (porto) 
P4 Via Orlando 307 608 P4 Via Orlando 230 138 
PS | ViaCateratte 344 98 PS | Via Cateratte 8 18 
P6 Via Mattei 263 412 P6 Via Mattei 82 80 
P7 Via Firenze 1790 846 P7 Via Firenze 167 97 
(sud) (sud) 

P8 | Via Enriques 45 255 P8 | Via Enriques SZ 75 
P9 Via Nord 246 484 P9 Via Nord 39 61 
P10 Via Pian di 561 632 P10 Via Pian di 56 121 

Rota Rota 
PII Connessione 793 756 PII Connessione 61 4 
Aurelia Aurelia 
pio | VIStuenZe: | 1185 | 1499 ppi SIEPneze ig 93 
(nord) (nord) 
TOTALE 7044 7044 TOTALE 1328 1328 
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Per poter ridistribuire questi flussi nella nuova rete è stato necessario fare 


delle ipotesi: 


e Peri veicoli leggeri 


. Come conseguenza del prolungamento di via 


Enriques il suo centroide sparisce per cui tutto il 
flusso dei veicoli leggeri entrante e uscente da 
questo è ridistribuito negli altri centroidi, in 


modo percentuale in base ai loro flussi. 


. Come conseguenza della divisione del traffico 


industriale da quello urbano nel nodo G, l'85% 
dei veicoli leggeri uscenti ed entranti da via 
Orlando sono stati ridistribuiti, in via Mattei e in 
via delle Cateratte, in modo percentuale in base ai 


loro flussi. 


. Diventando via delle Cateratte ad unico senso di 


marcia, entrante nella rete, il flusso uscente da 
questa è stato aggiunto al flusso uscente da via 


Mattei. 


. Avendo creato una nuova via per Tirrenia l'85% 


del traffico urbano che prima uscivano dal nodo 
FI-PI-LI (ovest) per andare in questa direzione è 
stato spostato nel nuovo nodo "Nuova via per 


Tirrenia". 
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e Peri veicoli pesanti 

1. Come conseguenza del prolungamento di via 
Enriques il suo centroide sparisce per cui tutto il 
flusso dei veicoli pesanti entrante e uscente da 
questo è ridistribuito negli altri centroidi, in 
modo percentuale in base ai loro flussi. 

2. Come conseguenza della divisione del traffico 
industriale da quello urbano nel nodo G, tutti i 
veicoli pesanti uscenti ed entranti da via Mattei e 
Via delle Cateratte sono stati aggiunti al flusso 


uscente ed entrante in via Orlando. 


Di seguito si riportano i valori dei flussi di traffico entranti e uscenti della rete 


di progetto: 
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Flussi di traffico entranti e uscenti 


della rete di progetto dei veicoli leggeri 


Flussi di traffico entranti e uscenti 
della rete di progetto dei veicoli pesanti 


(uvp/ora) (veic/ora) 
Centroide Entranti | Uscenti Centroide Entranti | Uscenti 
P1 | FI-PI-LI (Est) 867 826 P1 | FI-PI-LI (Est) 209 241 
FI-PI-LI FI-PI-LI 
P2 (Ovest) 72 78 P2 (Ovest) 178 190 
P3 Via L. da Vinci 170 165 P3 Via L. da Vinci 148 179 
(porto) (porto) 
P4 Via Orlando 48 95 P4 Via Orlando 333 250 
PS | ViaCateratte 494 0 PS | Via Cateratte 0 0 
P6 Via Mattei 378 1065 P6 Via Mattei 0 0 
P7 Via Firenze 1801 878 P7 Via Firenze 174 9g 
(sud) (sud) 
pg.| SOA via per euro ll duo pg SIONI ReR Lg 0 
Tirrenia Tirrenia 
P9 Via Nord 248 502 P9 Via Nord 41 65 
P10 Via Pian di 565 656 P10 Via Pian di 58 128 
Rota Rota 
PII Connessione 798 384 PII Connes sione 63 7g 
Aurelia Aurelia 
pig | Sete 1193 | 1555 pp attore 124 99 
(nord) (nord) 
TOTALE 7044 7044 TOTALE 1328 1328 
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3.1.3. Determinazione matrice 0/D 


Le matrici 0/D del nuovo scenario sono state ricavate con lo stesso 


procedimento visto nel capito 2.2.4.. Si riportano di seguito le iterazione e le 


matrici finali dei veicoli leggeri e dei veicoli pesanti calcolate separatamente. 


e Matrice O/D veicoli leggeri 


"Iterazione 0" 


* colonna + 


Au, 


+ 


826| 78 |165| 95 1065 | 878 | 440 | 502 | 656 | 784 | 1555 
0|0|)0|0 0 0|0|0|0|0 0 
Test di Arresto: 0,14 <0,01 NON VERIFICA 


X Riga Ae,i 
Ù Ù 
0 875 8 
0 81 9 
0 190 | 20 
0 54 6 
0 564 | 70 
0 371 -7 
0 1756 | -45 
0 441 31 
0 265 17 
0 588 23 
0 812 14 
0 1046 |-147 
0 
0 
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"Iterazione 1" 


X Riga Ae,i 

Ù Ù 

10 0 867 0 

0 0 TR 0 

2 0 170 0 

0 0 48 0 

5 0 494 0 

4 0 378 0 

21 0 1801 | 0 

4 0 410 0 

3 0 248 0 

6 0 565 0 

9 0 798 0 

15 0 1193 | 0 
» colonna + |832| 79 | 166| 96 | 0 |1070|877 |445 |506 | 662 | 791 | 1521 
Auj >| 6{1|1|1{0]|5]|-1|{5]|4]|6)|7 |-34 
Test di Arresto: 0,02 <0,01 NON VERIFICA 


"Iterazione 2" 


X Riga Ae,i 

I I 

0 867 0 

0 12 0 

0 169 -1 

0 48 0 

0 493 -1 

0 CIO -1 

0 1802 1 

0 409 -1 

0 249 1 

0 566 1 

0 799 I 

0 1193 0 
X colonna +-|826| 79 |166| 96 | 0 |1066]|878 | 439 |502 | 655 | 784 |1553 
Auj —>| 0 | 1 1 10 1 O0|-1|0]|-1| 0 | -2 

Test di Arresto: 0,006 <0,01 VERIFICA 
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e Matrice O/D veicoli pesanti 


"Iterazione 0" 


Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


* colonna + 


Au, 


"Iterazione 1" 


X Riga Ae, 
i Ù 
00 0 198 |-11 
00 0 178 0 
00 0 150 2 
00 0 304 | -29 
00 00 0 0 
00 00 0 0 
00 0|9 188 | 14 
00 00 0 0 
00 00 46 3 
00 0 | 3 62 4 
00 0 | 3 69 6 
00 O | 6 134 | 10 
241|190|179|250| 0 | 0 | 98 | 0 | 65 [128] 78 | 99 
—>|0|0]|0|0|0|0|{0{|0]|0|0|0]|0 
Test di Arresto: 0,102 <0,01 NON VERIFICA 
X Riga Ae, 
I 
00 0 | 11 0 
00 0|9 0 
00 07 0 
00 O | 18 0 
00 00 0 0 
00 00 0 0 
00 08 13 0 
00 00 0 0 
00 00 3 0 
00 0 | 3 4 0 
00 0 | 3 5 0 
00 O | 6 0 0 


00 
N 
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X colonna + | 240 |191|180|244| 0 | 0 |100] 0 | 66 |129| 79 100 
Auj —>|-1[{1|1|-6|0{|0|2|0|1|{1|{1}1 


Test di Arresto: 0,026 <0,01 NON VERIFICA 


"Iterazione 2" 


X Riga Ae, 
Ps | P6 |P7 [Ps |P9| dl 
00 0 noi 
0|0 0 [LS 
00 e: 
Op 
0 09 04 60800 0 
09 ROS (460 0 
07 F0-4] L09058 13 
0{|0/0|0/0 0 
On |04 34 OAK O 
004043 4 
0g f0g/054 045 9 
0704094 KOAE6 0 
© colonna — | | 0° | 
Au,j -1 -1|0)|0 00 


Test di Arresto: 0,008 <0,01 VERIFICA 
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3.1.4. Flussi di traffico nel nuovo scenario 


A questo punto le matrici trovate nel paragrafo precedente vengono inserite in 
Aimsun®M che, utilizzando l’algoritmo di Frank e Wolfe, ha determinato la 
distribuzione dei flussi negli archi della rete in condizioni di equilibrio 
deterministico, per i soli veicoli leggeri e i soli veicoli pesanti. 


Di seguito sono riportati i flussi ricavati per 1 vari quadri del modello: 


Quadro A 


Quadro E 
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Flussi quadro A 


[v] a 
n29 


via da Vinci 
(porto) 


Auto/h 
Pesanti/h 
pe/h 


Auto/h J FIPILI 
Pesanti/h 
pe/h 


via da Vinci Y 


[/ 


/ rotatoria di collegamento al 


“via Orlando nuovo Ponte Genova 


a 
Nuova via per Tirrenia e sii 
N39 Prali 


FIPILI Est 


Nuova via 
per Livorno 


nuovo collegamento 
con via Firenze 


via Enriques 


Flussi quadro B 
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a via Enriques 


Auto/h 
Pesanti/h 
pe/h 


via Firenze 47 
74 


nuova via 
per Tirrenia Auto/h 


e FIPILI Pesanti/h 
pe/h 


via Enriques 


nuova via” N 


per Tirrenia 
e FIPILI 


via Pian di Rota 


via Nord 
Flussi quadro D 
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Flussi quadro E 


\2/ Connessione con 
via Aurelia, A12, FIPILI 


Auto/h 
Pesanti/h 
pe/h 


3.2. Esame dettagliato di alcuni nodi 


Nel presente lavoro è stato scelto di esaminare le prestazioni dei due nodi più 
significativi: il nodo G e il nodo B. A questi è stato applicato il flusso della 
rete di progetto e sono state valutate le differenze prestazionali nel caso in cui 
si lasci l'intersezione come è attualmente oppure come è stata ipotizzata nella 


nuova rete. 
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3.2.1. Nodo G 


Per calcolare le prestazioni del nodo G sono stati utilizzati i flussi ricavati 
dalla simulazione di Aimsun®M per la nuova rete, ripresi dal quadro C del 


paragrafo 3.1.4. e riportati in figura 8 


96 auto/h 
249 pesanti/h 
594 pc/h 


Via L. da Vinci 
Via Orlando 


48 auto/h 
332 pesanti/h 
712 pc/h 


256 auto/h 
192 pesanti/h 
640 pc/h 


212 auto/h 
260 pesanti/h 


722 pc/h 
7 auto/h 
194 pesanti/h 
395 pc/h 
81 auto/h 286 auto/h 
108 pesanti/h 174 pesanti/h 


297 pc/ih 633 pc/h 


Nuovo ponte 
282 auto/h Genova 
189 pesanti/h 


660 pc/h 


870 auto/h 
0 pesanti/h 
870 pc/h 


Via Genova 


Via Mattei 


493 auto/h 
0 pesanti/h 
493 pc/h 


1066 auto/h 
0 pesanti/h 
1066 pc/h 


Via:delle Galeratio Rotatoria via Orlando 


Figura 8: flussi dei veicoli leggeri, pesanti e complessivi del nodo G della rete di progetto 


3.2.1.1. Rotatoria attuale 
Avendo calcolato il carico dalla simulazione di Aimsun®M usando lo schema 


della rete di progetto, per poterlo applicare alla rotatoria attualmente in uso, 


riportata in figura 9, è stato necessario formulare alcune ipotesi: 
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Figura 9: rotatoria attuale del nodo G 


e veicoli leggeri 
1. Non essendo la rotatoria attuale provvista di un'uscita 
per via Orlando il suo flusso uscente è stato sommato a 
quello uscente di via Mattei, che più a sud è collegata a 


via Orlando. 


94 


Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


e veicoli pesanti: 


di 


Poiché la rotatoria attuale non ha un collegamento 
diretto tra via Mattei, via delle Cateratte e il ponte tutto 
il loro flusso necessita di passare dall'approccio C. Per 
cui 1 flussi entranti e uscenti da via Mattei e via delle 
Cateratte sono stati aggiunti rispettivamente al flusso 


uscente e entrante di via Genova. 


. Essendo la rotatoria attuale sprovvista del bypass tra 


via Genova e via L. da Vinci, il flusso del bypass è 
stato aggiunto sia al flusso entrante di via Genova che 


al flusso uscente di via L. da Vinci. 


. Non essendo la rotatoria attuale provvista di un'uscita 


per via Orlando il suo flusso uscente è stato sommato a 
quello uscente di via Mattei, che più a sud è collegata a 


via Orlando. 


. Essendo la rotatoria attuale sprovvista del bypass tra 


via Genova e via L. da Vinci, il flusso del bypass è 
stato aggiunto sia al flusso entrante di via Genova che 


al flusso uscente di via L. da Vinci. 


Con queste considerazioni, dai flussi riportati in figura $8, sono stati ricavati i 


flussi entranti e uscenti espressi in uvp (1 veicolo leggero = 1 uvp; per 


uniformare i veicoli pesanti, si è utilizzato 1 veicolo pesante= 2 uvp) della 


rotatoria attuale che sono riportati nella tabella seguente: 
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Portate di traffico totali uvp/h 


Ramo Entranti |Uscenti 
Via Mattei SI 1660 
Via delle Cateratte 493 0 
Via Genova 1700 1530 
Via L. da Vinci 640 1592: 
Via Orlando 712 0 
TOTALE 3922 3922 


Da questi, come visto nel paragrafo 2.2.4., è stata ricavata la matrice O/D 


della rotatoria, riportata di seguito: 


Somma colonna 
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3.2.1.1.1. Metodo Inglese 


L'attuale rotatoria nel nodo G è stata costruita con il concetto delle rotatorie di 


1° generazione, ovvero dare la precedenza ai rami, e adattata a essere una 


rotatoria di 2° generazione con precedenza all'anello. 


Per questo motivo si utilizza, per il calcolo delle prestazioni, il metodo inglese 


che calcola 1l valore della massima portata oraria Qmax che può essere ospitata 


in un tratto di scambio e che è funzione della larghezza media all'inizio del 


tratto di scambio, della lunghezza del tratto e dell'incidenza proporzionale alle 


correnti di intreccio sul flusso di traffico. Il calcolo segue la formula: 


con: 


A = costante pari a 354 

w = larghezza della sezione stradale 
nel tratto di intreccio (m) 

L = lunghezza del tratto di intreccio 
(m) 

m= (mj + m-)/2 larghezza media 
nella sezione di imbocco del 


tratto di intreccio (m) 


b+c 


P= bici Proporzione delle 


correnti di scambio rispetto al 


flusso totale nel tratto di intreccio 


Channelization 
island 


Proportion of traffic weaving N 
b+c 
a+tb+c+d 


pa 


Note: This is the mirror image of the actual 
British left-hand driving situation. 
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Per calcolare i livelli di servizio dei tronchi di scambio si fa riferimento alla 
seguente tabella: 


<0,20 
0,20 + 0,50 


0,50 + 0,70 

0,70 + 0,85 

0,85 = 0,99 
1 


con Q;= a+b+c+d effettiva portata transitante. 


3.2.1.1.2 Prestazioni dell'intersezione 


Con i valori geometrici ripresi dalla figura 9, la matrice dei flussi ricavata nel 
paragrafo 3.2.1.1. e con le formule viste nel paragrafo precedente sono stati 
calcolati i livelli di servizio dei singoli tronchi di scambio e dell'intera 


intersezione: 


Tronco 1 (Via Orlando-Via Mattei) 


A | w|L|mlm2jmla |b c d | P|Qmax | Qt |Qt/Qmax| LOS 
354 |12,5|143|6,1] 9 |7,6|184|530|1476 |985/0,63|5152,50|3175| 0,62 


Tronco 2 (Via Cateratte-Via Genova) 


A | w|L|mlmM2|m})a)T]b c d P | Qmax | Qt |Qt/Qmax| LOS 
354 |11,5|128|5,7| 12 |8,9]309|185|1220]|1220|0,48|5554,92|2934| 0,53 
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Tronco 3 (Via Mattei-Via Cateratte) 


A 


ml 


m2 


mia b c d | P| Qmax 


Qt 


Qt/Qmax 


LOS 


302 


17 


8,3 


LI 


7,9] 0 | 379 | 0 |985 0,28 |3081,50 


1364 


0,44 


A 


ml 


m2 


mia b c d | P)| Qmax 


Qt 


Qt/Qmax 


302 


I,9 


(ie 


8,3 | 530 | 1161 | 202 | 127 0,67 | 2878,77 


2020 


0,70 


LOS dell'intera intersezione 


Qt/Qmax LOS 


0,60 


LOS 


99 


Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


3.2.1.2. Rotatoria di progetto 


Figura 10: rotatoria di progetto del nodo G 


Per calcolare le prestazioni della nuova rotatoria è stato usato il metodo HCM 


2016. 


3.2.1.2.1. HCM 2016 


Nell’Highway Capacity Manual 2016 viene proposto un modello analitico 
basato sulla teoria dell’accettazione del varco (Gap Acceptance). Dal 
momento che la rotatoria impone unicamente svolte a destra ai conducenti che 
dalle entrate si immettono nell’anello, le caratteristiche di accettazione del 
varco temporale sono assimilabili a quelle che si trovano nelle corsie di svolta 


a destra nelle intersezioni a priorità non semaforizzate. L’entità del varco è 
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strettamente dipendente dalla distribuzione degli spazi interveicolari nel 


flusso principale, nel caso delle rotatorie quello dell’anello. 
Due parametri definiscono la gap acceptance: 


t., l’intervallo critico. È il varco di tempo minimo perché un veicolo in 
attesa in un ramo di ingresso possa immettersi nell’anello. 

- t, il tempo di follow-up. È il tempo necessario affinché un secondo 
veicolo accodato al primo possa immettersi nella corrente veicolare 


dell’anello nello stesso varco del primo. 


Ne deriva che il tempo minimo necessario a due veicoli accodati per 
immettersi è dato dalla somma di intervallo critico e tempo di follow up 
(t.+t). Se 1 veicoli accodati sono più di due il tempo minimo diviene t:-+Nt; 
con N pari al numero di veicoli. 

Nei modelli più semplici la distribuzione degli arrivi dei veicoli circolanti 
nell’anello davanti ad un ramo d’entrata è di tipo normale, e 
conseguentemente la distribuzione dei varchi nel flusso è di tipo esponenziale. 
La probabilità che un varco sia minore di un determinato valore t è data 


dall’equazione: 


PIT<t]=e-® 


in cui 


A = valore medio degli arrivi per il flusso di disturbo [veic/h] 


t = dimensione headway [h] 
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Combinando la formula esposta con t; e tr, si arriva alla formula HCM per la 
stima della capacità di un’entrata (formula per una rotatoria con corsia singola 


nell’anello). 


-0ce'tc 
Qo . el 3600 ) 
Co = TL *t 
Qete 
1 —_ el 3600 ) 


dove: 

C.= capacità dell'approccio [veic/h] 
Q; = flusso circolante [veic/h] 

t;= gap critico [s] 


t: = tempo di follow up [Ss] 


Il modello generale deve essere opportunamente calibrato scegliendo i valori 
di t, e di t; opportuni in base al comportamento degli utenti e alla geometria 
della rotatoria. Secondo le osservazioni dei ricercatori statunitensi sono stati 
individuati i valori di massimo e di minimo per l’intervallo critico e il tempo 


di follow up, valori riportati nella tabella sottostante. 


Massimo (s) Minimo (s) 
Intervallo critico t; 4.6 4.1 
Tempo di follow up tr 3.1 2.6 


Sostituendo i valori nell’Equazione precedente si ottiene il modello 
semplificato per la capacità di un’entrata per una rotatoria ad una corsia 


nell'anello: 


C, = 1380 + e(-102:1079)0c 


Per rotatorie con due corsie nell’anello vengono proposte le seguenti formule: 
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“entrata a corsia singola C, = 1420 + e(-0.85:107)0c 


-entrata a corsia doppia Cop = 1420 - e(-0.85:107°)Qc 


Co, = 1380: e(-0.92:107)Qc 
Per le corsie di bypass viene proposta le seguente formula: 


Cebypass = 1420: e(-0.85:107?)Qc 


Nel Paragrafo 4 del Capitolo 22 dedicato alle rotatorie del HCM 2016, si 
introduce la possibilità di poter generalizzare il modello di capacità usando le 
seguenti formule: 

CsA 0 
In cui: 


DA 
2330001 (e a) 
tr “3600 


Con il modello generalizzato è possibile inserire l’intervallo critico t; e il 
tempo di follow up t; idonei in modo che 1 risultati ottenuti siano aderenti 
all’effettivo comportamento degli utenti a livello locale. I valori ottenuti da 
varie campagne di osservazione a livello locale hanno portato alla 
determinazione di valori per l’intervallo critico ed il tempo di follow up che si 
discostano, talvolta in maniera marcata, da quelli determinati dallo studio 


HCM. I dati ricavati possono essere confrontati nella tabella riportata. 
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_____  b@l@la[ B_ | Gem | 
Single Lane Roundabout 

2,50] 4.90] 1440 

2,6 1364*exp((-0,70*104-3)#Q.) 
Left Lane Roundabout 

1130*exp((-0,75*109-3)*Q%) 

NCHRP 572 4,50 

1390*exp((-0.71*101-3)*Q%) 
Right Lane Roundabout 

[NCHRP 572___|3,10|4,20|1161|0,000736] 1161*exp((-0.736*104-3)*Q.)_| 


Nell'HCM la metodologia usata per il calcolo della capacità di un approccio e 
la sua valutazione in termini di livello di servizio si articola in dodici fasi, 
schematicamente riassunte nella figura seguente, e mostra come procedere 
dalla determinazione dei flussi fino all’inquadramento del livello di servizio 


di ogni singolo approccio. 
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Step 1: Convert movement demand volumes to flow rates 


Step 2: Adjust flow rates for heavy vehicles 


Step 3: Determine circulating and exiting flow rates 


Step 4: Determine entry flow rates by lane 


Step 5: Determine the capacity of each entry lane and 
bypass lane as appropriate in passenger car equivalents 


Step 6: Determine pedestrian impedance to vehicles 


Step 7: Convert lane flow rates and capacities into 
vehicles per hour 


Step 8: Compute the volume-to-capacity ratio for each lane 


Step 9: Compute the average control delay for each lane 


Step 10: Determine LOS for each lane on each approach 


Step 11: Compute the average control delay and determine 
LOS for each approach and the roundabout as a whole 


Step 12: Compute 95th percentile queues for each lane 


Il livello di servizio viene valutato in base al tempo che un veicolo deve 
attendere prima di potersi immettere nell’anello, secondo gli intervalli 


riassunti nella figura seguente. 
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Tempo di attesa (s) | LOS | Colore | Deflusso 
0-10 A Ottimo 
> 10-15 B Discreto 
> 15-25 C Buono 
> 25-35 D Sufficiente 
> 35-50 E Insufficiente 
> 50 F Congestione 


Una volta calcolata la capacità con il metodo appropriato per il tipo di 
rotatoria, per stimare i tempi di attesa e la lunghezza del 95° percentile 
vengono impiegate le formule proposte dall’HCM. Il modello per la stima del 
tempo di attesa è stato desunto considerando la rotatoria come un’intersezione 
a priorità, ipotizzando una distribuzione di Poisson per gli arrivi e una 


disciplina della coda di tipo FIFO (First In First Out). 


Il tempo di attesa su una singola corsia di entrata viene stimato con la 


seguente equazione: 


(0, 
4507 


3600 


d=—_+9007 x-14+ [(x-1)?+ + 5 min[x, 1] 


In cui: 


d = tempo di attesa medio dei veicoli nell’approccio [s/veic] 


c = capacità [veic/h] 


x = indice di saturazione == 


Qe (Portata entrante 
Ce 


) l'HCM prevede, per una 


capacità 


rotatoria a due corsie in ingresso, in mancanza 
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di dati di suddividere il flusso entrante (Qe) in 
53% per la corsia di sinistra e 47% per quella 
destra. 


T = periodo di riferimento [h] (T è pari a 1 per analisi di un’ora; 0.25 per 


analisi di un quarto d’ora) 


Una volta stimati i tempi di attesa su entrambe le corsie di un approccio si 
calcola il tempo di attesa medio dell’approccio stesso eseguendo una media 


pesata dei tempi delle singole corsie, usando come pesi 1 flussi entranti. 


Asx* Usx + dax' Vax + diypassVbypass 
dapproccio 7 


Usx st Vax 2 Ubyapass 
In cui 

d,, = tempo di attesa sulla corsia sinistra 

dix = tempo di attesa sulla corsia destra 

diypass = tempo di attesa sul bypass 

Vv, = flusso entrante sulla corsia sinistra 

Vax = flusso entrante sulla corsia destra 


Vbypass = flusso entrante sul bypass 


L’equazione che permette di ricavare il valore del tempo di attesa medio 
dell’intera rotatoria è una media pesata del tempo di attesa dei vari approcci. I 
pedici si riferiscono ai rami in ingresso: 


d _ Xd;v; 
tatoria T 
rotatoria Y v; 
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Il 95° percentile della coda per una corsia viene stimato attraverso 
l’equazione: 


(0, 
Qos = 9007 |x-1+ |(x- 1)2 + SUVA t (=) 


In cui: 
Q95 = 95° percentile della coda [veic] 


c = capacità [veic/h] 


Sa ; , portata entrante 
x = Indice di saturazione (—__ 


) 


capacità 


T = periodo di riferimento [h] (T è pari a 1 per analisi di un’ora; 0.25 per 
analisi di un quarto d’ora) 


Se si desidera esprimere la lunghezza della coda in metri anziché in veicoli è 
sufficiente moltiplicare il valore ottenuto per 6 metri (vale a dire la lunghezza 
media di un veicolo, considerando anche la distanza di rispetto dal veicolo 


che segue) e si divide per il numero di corsie considerate n: 


Q95 ‘6 
L = 
q95 n [m] 
3.2.1.2.2. Prestazione dell'intersezione 


Con i flussi di figura 8, come visto nel paragrafo 2.2.4., è stata ricavata la 


matrice O/D della rotatoria riportata di seguito. 
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Somma colonna na 594 660 395 1649 


Con le formule dell'HCM 2016 viste nel paragrafo precedente sono stati 


calcolati i livelli di servizio dei singoli approcci e dell'intera intersezione: 


Approccio 1 via Orlando 


Delay 
. Delay ; LOS 
Corsia | Ve Ce Qe x |T (sec) Qos | Los e sppricdia 
Sx |416|1034,53| 378,42 10,37| 1 7,31 2 12 701 
Dx |416|1046,39| 335,58 |0,32| 1 6,67 2 12 ” 
Approccio 2 via Genova 
Delay 
è Delay è LOS 
Corsia | Vc Ce Qe x |T (660) Qos | Los 0 daprocaa 
Sx |339|1092,66| 156,88 |0,14| 1 4,56 1 6 
Dx |339|1099,75| 139,12 |0,13| 1 4,38 1 6 319 
Bypass | 395 | 1015,02 | 336,00 | 0,33 | 1 6,95 2 12 
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Approccio 3 via L. da Vinci 


Delay 
i Delay î LOS 
Corsia | Ve Ce Qe x T (#ec) Q95 | Los ara sopra 
Sx |240|1172,23| 338,67 |0,29| 1 5,76 2 12 5.60 
Dx |240|1172,39| 300,33 | 0,26 | 1 5,41 2 12 i 


Tempo di attesa 
nell'intersezione (s) 


LOS Intersezione 


6,17 


3.2.1.3. Confronto tra le prestazioni delle rotatorie 


Di seguito si riportano a confronto le prestazioni delle due rotatorie analizzate 


Rotatoria Attuale Nuova rotatoria 
Qt/Qmax LoS Tempo di attesa (s) LoS 
0,60 6,17 


Dal risultato si può notare come la rotatoria di progetto abbia un livello di 


servizio superiore rispetto a quella attuale a parità di flusso. 
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3.2.2. Nodo B 


Per calcolare le prestazioni del nodo G sono stati utilizzati i flussi ricavati 


dalla simulazione di Aimsun®M della nuova rete, riportati in figura 11. 


1460 auto/h 
85 pesanti/h 
1630 pc/h 


Rotatoria Via Firenze 
632 pcih AREA 


Figura 11: flussi dei veicoli leggeri, pesanti e complessivi del nodo B della rete di progetto 
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3.2.2.1. Intersezione semaforica 


Avendo calcolato il carico dalla simulazione di Aimsun®M usando lo schema 
della rete di progetto, per poterlo applicare all'intersezione semaforica 
attualmente in uso, riportata in figura 12, è stato necessario formulare alcune 


ipotesi: 


a 


Figura 12: attuale intersezione semaforica del nodo B 
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e veicoli leggeri 


e veicoli pesanti: 


. Essendo l'attuale intersezione semaforica sprovvista del 


bypass tra via Enriques e via Firenze, il flusso del 
bypass è stato aggiunto sia al flusso entrante di 


Enriques che al flusso uscente di via Firenze. 


. Essendo l'attuale intersezione semaforica sprovvista del 


bypass tra via Firenze e via Enriques, il flusso del 
bypass è stato aggiunto sia al flusso entrante di via 


Firenze che al flusso uscente di via Enriques. 


. Essendo l'attuale intersezione semaforica sprovvista del 


bypass tra via Enriques e via Firenze, il flusso del 


bypass è stato aggiunto sia al flusso entrante di 


Enriques che al flusso uscente di via Firenze. 


. Essendo l'attuale intersezione semaforica sprovvista del 


bypass tra via Firenze e via Enriques, il flusso del 
bypass è stato aggiunto sia al flusso entrante di via 


Firenze che al flusso uscente di via Enriques. 


Con queste considerazioni, dai flussi riportati in figura 11, sono stati ricavati i 


flussi entranti e uscenti espressi in uvp (1 veicolo leggero = 1 uvp; per 


uniformare i veicoli pesanti, si è utilizzato 1 veicolo pesante= 2 uvp) 


dell'attuale intersezione semaforica che sono riportati nella tabella seguente. 
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Portate di traffico totali uvp/h 
Ramo Entranti |Uscenti 
PI Via Firenze Ovest 2377 1710 
P2 Via Nord 331 421 
P3 Via Firenze Est 1203 1632 
P4 Via Enriques 1171 1319 
TOTALE 5082 5082 


Da questi, come visto nel paragrafo 2.2.4., è stata ricavata la matrice O/D 


dell'intersezione semaforica riportata di seguito. 


Somma colonna - 1710 421 1631 1319 5081 
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3.2.2.1.1. HCS+ 


Per l'intersezione semaforica è stato adottato il software HCS+ che calcola la 
capacità e i livelli di servizio delle infrastrutture di trasporto come da Manuale 
HCM. 
Per poter usare il programma HCS+ per le intersezioni semaforiche è 
necessario conoscere: 
1. La geometria ed 1 flussi: i primi sono stati ricavati durante la campagna 
di rilievo, mentre i secondi sono stati ottenuti 
dalla matrice nel paragrafo precedente. 
2. Il ciclo semaforico e le manovre: Sono stati ricavati durante la 
campagna di rilievo. A seguito 
riportiamo un'immagine con i 


tempi del semaforo e le manovre. 


N di 
J i 
j è RK È =A VA 


Fase A Fase B Fase C 
Tempo di verde 30 (s)| Tempo di verde 15 (s)| Tempo di verde 15 (s) 
Tempo di giallo 5 (s) | Tempo di giallo 5 (s) | Tempo di giallo 5 (s) 


Figura 13: ciclo e manovre del semaforo 
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I dati vengono inseriti nel programma in apposite celle, di seguito si riportano 


alcune immagini della procedura: 


e Inserimento dei dati personali e dell'intersezione: 


Date [26/10/2017 


Analyst [DavidZanobeti 00° 
| Agency/Co IUMIPPO 
Units: U. S. Metric 
Analysis Time Period [AM peak 

Project Description [Tesi di Laurea Magistrale 


| => - : 
East/Aw'est Street Name [via Firenze North#South Street Name (via Enriquez 


Intersection 
Area Type 
Jurisdiction 


Analysis Year 


{via Firenze - Via Enriquez 


TT CBD or Similar 


[Livoro 
fa01 7? 


e Inserimento della geometria dell'intersezione 


Sl 
EI 
EJ 
dl 
si 
Al 
al 


Ai] San] 8] Da] na] na) pe 


isso 


Revert 


Cancel 


[n ti ls n lin fe be 
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e Inserimento dei flussi per ciascuna corsia 


-GEOMETRY and VOLUME È È: 
| Eastbound | Westbound | Northbound | Southbound 


| Left | Thru | Right | Left | Thru | Right Left | Thru | Right | Left | Thru | Right 


Number of Lanes and Usage 
fio efifo eat ai ee St a e o 
Shared | Shared | Shared | Shared Shared | Shared Shared | Shared 


JE T [È TR L TR L TR 
Receiving Lanes Ali Als fs hi. “ 
Volume (ph), Increment All fo % 4|- | Duration  |0.25 hours 


ross fina fo ana Ros 7 fo fe po fre 
Peak Hour Factor, PHF, All [0.30 — 
[OS TS Ts fo oe oe iS Te osp po i 
Peak-15 Minute Volume [+] 

263 329 0 55 195 54 52 15 18 62 51 182 


e Inserimento delle fasi semaforiche e dei tempi di verde e di giallo 


PHASING DESIGN 
PRESET PHASING 


EW' [Select New Preset “| 
NS [Select New Preset <] 


[Select Actuation Type v ] 


Cycle Length [75.0 Estimation/Optimization | 
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Una volta inseriti i dati è possibile stampare un report nel quale si trovano vari 
valori tra cui: la capacità delle corsie (C), l'indice di saturazione (x) ed inoltre 
il tempo di attesa (sec) e il livello di servizio per ogni corsia, per ogni 
approccio e per l'intera intersezione. 


L'immagine seguente mostra la parte del report che fornisce 1 dati ricavati: 


Intersection Performance Summary 
Appr/ Lane Adj Sat Ratios Lane Group Approach 
Lane Group Flow Rate 
Grp Capacity (s) v/c g/C Delay LOS Delay LOS 


Eastbound 

L 357 1785 2.94 0.20 909.8 F 

LD 1431 3577 0.92 0.40 322 € 421.8 F 
Westbound 

L 357 1785 0.61 0.20 34.9 C 

TR 1385 3462 0.72 0.40 22-% € 24.4 € 
Northbound 

L 100 501 2.05 0.20 535.4 F 

TR 345 1724 0.37 0.20 28.9 € 341.6 F 
Southbound 

L 254 1269 0.98 0.20 80.6 FF 

TR 332 1659 2.80 0.20 847.5 F 685.9 F 


Intersection Delay = 383.1 (sec/veh) Intersection LOS = F 


HCS+: Signalized Intersections Release 5.1 


3.2.2.1.2. Prestazioni dell'intersezione 


Inseriti i flussi ricavati dalla matrice del paragrafo 3.2.2.1., 1 valori 
geometrici, il ciclo semaforico e le manovre mostrati in figura 13 nel 
programma HCS+ sono stati calcolati i livelli di servizio di cui riportiamo in 
seguito delle tabelle riassuntive. 
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Approccio 1 via Firenze Ovest (Estbound) 


Corsia V C x Delay Delay LOS i 
(sec) approccio (s) | approccio 
Left 1049 | 357 |2,94| 909,80 
Thu |1313|1431|0,92| 32,20 iii 
corsie) 
Approccio 2 via Nord (Northbound) 
Corsia V C x 200) DERE LS » 
(sec) approccio (s) | approccio 
Left 205 | 100 |2,05| 535,40 341.65 
Thru/Right | 127 | 345 |0,37| 28,90 i 


Approccio 3 via Firenze Est (Westbound) 


Corsia V C x pene delay LOS . 
(sec) approccio (s) | approccio 
Left 218 | 357 0,61) 34,90 
i 24,41 
Thn/Right | 992 |1385|0,72] 22,10 i 
(2 corsie) 
Approccio 4 via Enriques (Southbound) 
Corsia V EC x DAY DEAN LOS i 
(sec) approccio (s) | approccio 
Left 248 | 254 |0,98| 80,60 685.91 
Thru/Right | 929 | 332 |2,80| 847,50 i 


Tempo di attesa 
nell'intersezione 


(8) 


383,18 


LOS Intersezione 
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Dai risultati precedenti si vede come l'intersezione risulti congestionata. Per 
valutarne il funzionamento sono state applicate varie condizioni di flusso, 1l 
90%, l'80% e il 70% del totale. Di seguito si riportano le tabelle riassuntive 


dei vari casi di carico: 


e Casoconil flusso pari al 90% del totale 


Approccio 1 via Firenze Ovest (Estbound) 


Delay Delay LOS 
(sec) approccio (s) | approccio 
Left 944 | 357 |2,64| 777,90 


Corsia V C x 


Thru(@ |1182|1431|0,83| 25,70 “ù 
corsie) 
Approccio 2 via Nord (Northbound) 
Corsia V C x DElaY Delay LO6 i 
(sec) approccio (s) | approccio 
Li 185 | 100 |1,85| 448,30 287.64 
Thru/Right | 114 | 345 [0,33 28,30 


Approccio 3 via Firenze Est (Westbound) 


Corsia v c > Delay Delay LOS 
(sec) approccio (s) | approccio 
Left 196 | 357 |0,55| 32,90 
i 24,54 
Thru/Right | g93 |1385|0,64| 22,70 ? 
(2 corsie) 
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Approccio 4 via Enriques (Southbound) 


Corsia V C x DS Door LO5 . 
(sec) approccio (s) | approccio 
Left 223 | 257 |0,87| 59,40 582.39 
Thru/Right | 836 | 332 |2,52| 722,00 ” 
Tempo di attesa 
nell'intersezione LOS Intersezione 
(5) 
326,79 
e Caso conil flusso pari al 80% del totale 
Approccio 1 via Firenze Ovest (Estbound) 
Corsia V E x coni Delay 05 î 
(sec) approccio (s) | approccio 
Left 839 | 357 |2,35| 646,20 
299,4 
Thu | 1050 |1431|0,73| 22,50 si 
corsie) 
Approccio 2 via Nord (Northbound) 
Corsia V C x DAY Delay LOS ; 
(sec) approccio (s) | approccio 
Left 164 | 100 |1,64| 358,40 231.90 
Thru/Right | 102 | 345 |0,29| 27,70 ” 
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Approccio 3 via Firenze Est (Westbound) 


Corsia V C x Delo) DEI Lo ; 
(sec) approccio (s) | approccio 
Left 174 | 357 |0,49| 31,30 
i 21 
Thrw/Right | 794 |1385/0,57| 19,20 ne 
(2 corsie) 


Approccio 4 via Enriques (Southbound) 


Corsia v c i Delay Delay LOS l 
(sec) approccio (s) | approccio 
Li 198 | 260 [0,76] 47,10 480,90 
Thru/Right | 743 | 332 |2,24| 596,70 


Tempo di attesa 
nell'intersezione LOS Intersezione 


(Ss) 
270,87 


e Casoconil flusso pari al 70% del totale 


Approccio 1 via Firenze Ovest (Estbound) 


Delay Delay LOS 
(sec) approccio (s) | approccio 
Left 734 | 357 |2,06| 514,80 


Thru (2 | 919 | 1431|0,64| 20,40 29221 
corsie) 


Corsia V C x 
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Approccio 2 via Nord (Northbound) 


Corsia V C x 2099 Doo Lo ; 
(sec) approccio (s) | approccio 
Left 144 | 100 |1,44| 275,50 180,48 
Thru/Right |88,9| 345 |0,26| 27,10 ” 
Approccio 3 via Firenze Est (Westbound) 
Corsia V C x Delay DEIAY LOS ì 
(sec) approccio (s) | approccio 
Left 153 | 357 [0,43 30,00 
ThruRight | 694 | 1385 |0,50| 18,20 sii 
(2 corsie) 
Approccio 4 via Enriques (Southbound) 
Corsia V C x I0Ry Day LOS ? 
(sec) approccio (s) | approccio 
174| 2 4 
Li 63 | 0,66 0,00 380,82 
Thru/Right | 650 | 332 |1,96| 471,80 


Tempo di attesa 
nell'intersezione 


(8) 


216,40 


LOS Intersezione 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


3.2.2.2. Rotatoria di progetto 


Figura 14: rotatoria di progetto del nodo B 


Per calcolare le prestazioni della nuova rotatoria è stato usato il metodo HCM 


2016 visto nel paragrafo 3.2.1.2.1. 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


3.2.2.2.1 Prestazioni dell'intersezione 


Con i flussi di figura 11, come visto nel paragrafo 2.2.4., è stata ricavata la 


matrice O/D della rotatoria riportata di seguito. 


Somma colonna 


Con le formule dell'HCM 2016 sono stati calcolati i livelli di servizio dei 


singoli approcci e dell'intera intersezione: 


Approccio 1 via Firenze Ovest (Estbound) 


n Qos | Los 
Sx 541 | 946,67 | 1259,81 | 1,33) 1 |619,14|168 | 1008 
Dx 541 |965,21 | 1117,19|1,16| 1 |317,33| 94 | 564 473,13 

Bypass | 852 | 688,29 851 1,24| 1 |\461,53| 95 | 570 


Delay LOS 


Corsia | Ve Ce Qe x |T È . 
approccio (s) | approccio 


Approccio 2 via Nord (Northbound) 


5 Del Del LOS 
SODO ce Qe = L. Qus | Los mi (s) | approccio 
Sx |2499 | 235,75 | 175,43 |0,74| 1 | 59,79 | 8 48 
Dx |2499|272,46 | 155,57 |0,57| 1 | 33,27 | 4 24 30,69 
Bypass | 419 | 994,52 221 0,22 1 | 5,76 1 6 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


Approccio 3 via Firenze Est (Westbound) 


è Del Del L 
Corsui.|; VG Ce Qe 3. (sec) Qos | Los eat (8) ni 
Sx 1193 | 595,87 | 481,24 |0,81| 1 | 33,82 | I1 | 66 
Dx |1193|633,44 | 426,76 |0,67| 1 | 20,56 | 6 36 24,39 
Bypass | 1028 | 592,65 295 |0,50| 1 | 14,54 | 3 18 


Approccio 4 via Enriques (Southbound) 


Delay 
(sec) Q95 | Los 


Sx _|1073|648,87| 169,6 |0,26| 1 | 8,82 | 2 12 
Dx |1073|684,49| 150,4 |0,22| 1 | 7,84 1 6 


Delay LOS 


Corsia | Ve Ce Qe x |T - . 
approccio (s) | approccio 


8,36 


Tempo di attesa 
nell'intersezione (s) 


LOS Intersezione 


308,82 


Dai risultati precedenti si vede come la rotatoria risulti congestionata. Per 
valutarne il funzionamento sono state applicate varie condizioni di flusso, il 
90%, l'80% e il 70% del totale. Di seguito si riportano le tabelle riassuntive 


dei vari casi di carico: 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


e Casoconil flusso pari al 90% del totale 


Approccio 1 via Firenze Ovest (Estbound) 


. Delay Delay LOS 
Coma e Se Qe i (sec) Qos | Los approccio (s) | approccio 
Sx 487 | 983,74 |1133,83 [1,15 308,60 | 93,3 | 560 
Dx |487 | 999,54 | 1005,47 | 1,01 94,87 |40,3| 242 190,62 
Bypass | 767 | 739,99 | 765,9 |1,04 141,67 | 41 | 246 
Approccio 2 via Nord (Northbound) 
è Delay Delay LOS 
Corsia Xe ce Qe z (sec) Qos | Los approccio (s) | approccio 
Sx |2249| 281,52 | 157,887 |0,56 31,60 |3,62|21,7 
Dx |2249]| 320,19 | 140,013 |0,44 22,09 |2,27|13,6 18,39 
Bypass | 377 | 1030,58 | 198,9 |0,19 5,29 |0,72| 4,3 
Approccio 3 via Firenze Est (Westbound) 
i Delay Delay LOS 
Cosa SS ce Qe 5 (sec) Qus | Los approccio (s) | approccio 
Sx |1074| 648,55 |433,116 |0,67 19,85 |5,73|34,4 
Dx |1074| 684,18 |384,084 |0,56 14,75 |3,75|22,5 15,99 
Bypass | 925 | 646,77 | 265,5 |0,41 11,48 |2,07|12,4 
Approccio 4 via Enriques (Southbound) 
. Delay Delay LOS 
Corn Ve ce Se È (sec) Qos | Los approccio (s) | approccio 
Sx 966 | 700,23 | 152,64 |0,22 7,66 |0,83| 5 739 
Dx |966 | 733,62 | 135,36 (0,18 6,94 |0,68 | 4,06 ” 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


Tempo di attesa 
nell'intersezione (s) 


126,52 


LOS Intersezione 


e Caso conil flusso pari al 80% del totale 


Approccio 1 via Firenze Ovest (Estbound) 


Delay 


Delay LOS 
(sec) Q9s Los 


Corsia | Ve Ce Qe x |T . a 
approccio (s) | approccio 


Sx | 433 | 1022,26 | 1007,85 |0,99| 1 | 75,83 |35,4| 213 


Dx | 433 | 1035,09 | 893,752 |0,86| 1 | 28,11 |15,5|93,3 48,11 


Bypass | 682 | 795,56 | 680,8 |0,86| 1 | 33,33 | 14,3 | 85,5 


Approccio 2 via Nord (Northbound) 


Delay 


Delay LOS 
(sec) Q9s Los 


Corsia | Ve Ce Qe x |T 2 . 
approccio (s) | approccio 


Sx |1999| 336,17 | 140,344 |0,42| 1 | 20,42 | 2,1 | 12,6 


Dx |1999| 376,29 | 124,456 |0,33| 1 | 15,93 |1,47|8,79 12,93 


Bypass | 335 | 1067,95 | 176,8 |0,17] 1 | 4,87 |0,59 3,57 


Approccio 3 via Firenze Est (Westbound) 


Delay 
(sec) 


Delay LOS 


Corsia | Ve Ce Qe x |T ; è 
approccio (s) | approccio 


Sx | 954 | 705,87 | 384,992 |0,55| 1 | 13,90 |3,52|2],1 


Dx | 954 | 739,00 | 341,408 |0,46| 1 | 11,35 |2,54|15,3 11,87 


Bypass | 822 | 705,83 236. |0,33|1 | 9,33 |.L5 |9,98 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


Approccio 4 via Enriques (Southbound) 


2 Delay Delay LOS 
po Ae ia Qe FI (sec) Qos | Los approccio (s) | approccio 
Sx 858 | 755,66 | 135,68 |0,18| 1 | 6,70 [0,66 3,93 6.45 
Dx | 858 | 786,28 | 120,32 [0,15 1| 6,17 |0,54|3,25 i 
Tempo ci de LOS Intersezione 
nell'intersezione (s) 
34,62 
e Casoconil flusso pari al 70% del totale 
Approccio 1 via Firenze Ovest (Estbound) 
ì Delay Delay LOS 
A pia Qe x | (sec) Qos | Los approccio (s) | approccio 


Sx 379 | 1062,29 | 881,867 |0,83| 1 | 23,29 |12,8| 77 


Dx | 379 | 107191 | 782,033 |0,73| 1 | 15,90 |7,69|46,1 19:13 


Bypass | 596 | 855,32 | 595,7 |0,70] 1 | 17,18 [6,55 | 39,3 


Approccio 2 via Nord (Northbound) 


. Dela Dela LOS 
Le De ce Qe = 4 (eo) Qos | Los sia (s) | approccio 
Sx |1749]| 401,44 | 122,801 |0,31| 1 | 14,44 |1,31|7,86 
Dx |1749| 442,21 | 108,899 |0,25| 1 | 12,03 |0,97| 5,85 STI 
Bypass | 293 | 1106,67 | 154,7 |0,14| 1| 4,48 |0,49|2,92 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


Approccio 3 via Firenze Est (Westbound) 


; Delay Delay LOS 
cora Ve Ce Qe = (sec) Qos | Los approccio (s) | approccio 
Sx 835 | 768,27 | 336,868 (0,44 10,52 |2,32|13,9 
Dx |835 | 798,20 | 298,732 [0,37 9,07 |1,78|10,7 9,32 
Bypass | 720 | 770,28 | 206,5 |0,27 7,72 |1,09| 6,57 
Approccio 4 via Enriques (Southbound) 
2 Delay Delay LOS 
Sora. We ti Qe è (sec) Qus | Los approccio (s) | approccio 
Sx 751 | 815,48 | 118,72 (0,15 5,89 |0,51|3,06 571 
Dx | 751 | 842,71 | 105,28 {0,12 5,51 |0,43|2,57 ” 


Tempo di attesa 
nell'intersezione (s) 


15,35 


LOS Intersezione 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


3.2.2.3. Confronto tra l'intersezione semaforica e la rotatoria 


Di seguito si riportano a confronto le prestazioni globali e di ogni singolo 


approccio dell'intersezione semaforica e della rotatoria analizzate: 


Confronto tra i singoli approcci del semaforo e della rotatoria 


0, 7 7 
%o W (Estbound) | S(Northbound) | E(Westbound) | N (Southbound) 


Los LoS Los 


341.65 24.41 6 
30,69 24,39 


Flusso 


Semaforo si 

= 
= Dr) 
Rotatoria | 19,13 9,32 [sm 


Confronto tra l'intersezione semaforica e la rotatoria 


cai 
Day © Day © 
(*] 


308.82 


Ss 
6.45 
380,82 
5.71 


. Ss 
2 18,39 
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Progetto di una nuova rete viaria a Nord-Ovest del porto di Livorno 


Dai risultati si può notare come la nuova rotatoria riesca a far defluire in 
modo migliore il traffico, infatti già con una riduzione del 20% del traffico 
totale ha un deflusso sufficiente. Mentre l'intersezione semaforica risulta 
sempre pesantemente congestionata anche con una riduzione del flusso pari al 


30% del totale. 
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4. Conclusioni 


Il quadro dell'attuale rete infrastrutturale della zona Nord-Ovest del porto di 
Livorno, ci restituisce una commistione tra il traffico urbano e il traffico 
industriale, come è stato possibile constatare dai rilievi e dalla simulazione 
effettuata sulla rete attuale. 

La simulazione sulla rete di progetto ha portato ad uno scenario di flusso 
diverso da quello della rete attuale, con una differente distribuzione del 
traffico industriale e del traffico urbano sugli archi della rete. Di seguito si 
riporta una tabella riassuntiva e gli istogrammi che mettono a confronto i 


flussi dei veicoli leggeri e pesanti degli archi della rete comuni tra i due 


scenari: 
Variazione di 
i ? Flusso flusso della rete 
Tipologia i 
Arco veicoli ; di progetto 
Rete Rete di rispetto alla rete 
Attuale Progetto attuale 


Svincolo FI- 
PI-LI - nodo 
FI-PI-LI 
(ovest) 


Leggeri 


FI-PI-LI 
Pesanti 


Nodo G - Leggeri 
Svincolo FI- 
PI-LI Pesanti 


Via L. Da 
Vinci 


Via Leggeri 
Firenze 
(Est) Pesanti 


Via Nodo B - Leggeri 
Firenze | Semaforo via 
(Ovest) Genova Pesanti 

Via 
Genova 

e via 
Mattei 


Semaforo via | Leggeri 103,36% 
Genova - via 
Mattei Pesanti -100,00% 
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Flusso veicoli leggeri 


M Flusso rete attuale 


M Flusso rete di Progetto 


Svincolo Nodo G- NodoA- NodoB- Semaforo 
FI-PI-LI - Svincolo  NodoB Semaforo via Genova 
nodo FI-PI-  FI-PI-LI via Genova — -via 

LI (ovest) Mattei 


Flusso veicoli pesanti 


Mm Flusso rete attuale 


M Flusso rete di Progetto 


FI-PI-LI Via L. Da Via Firenze Via Firenze Via Genova 
Vinci (Est) (Ovest) e via 
Mattei 
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Conclusioni 


Come si può notare dai risultati della tabella e dai grafici precedenti il flusso 
dei veicoli leggeri è diminuito sia in via L. da Vinci (-63,41%), sia nella parte 
terminale della FI-PI-LI (-84,56%), mentre il flusso dei veicoli pesanti è 
diminuito in Via Firenze (-43,84% e -10,86%), via Genova e via Mattei 
(-100%). 

I risultati ottenuti mettono in evidenza che via L. da Vinci diventa 
principalmente viabilità industriale infatti la maggior parte del traffico urbano 
si è spostato soprattutto in via Firenze e in via Enriques trasferendo il limite 
portuale da via L. da Vinci a via Enriques. Inoltre grazie al nuovo arco tra il 
nodo G e il nodo H si va a valorizzare la zona tra via L. da Vinci e via 
Enriques. 

Quindi gli obbiettivi prefissati per la rete di progetto esposti nel paragrafo 3.1 


sono stati confermati dai risultati ottenuti dalla simulazione. 
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